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ABSTRACT 

The  self-excited  oscillation  of  a  square  section  vertical 
cantilever  beam  and  a  square  section  sagging  suspended  cable  in 
a  uniform  wind  stream  was  investigated  theoretj  •'ally  and  experi¬ 
mentally  . 

Experimental  quasi -static  representation  of  aerodynamic 
forces  was  used  to  determine  the  responses  theoretically.  The 
excitation  is  due  to  negative  damping  type  aerodynamic  forces 
and  the  amplitude  is  limited  as  the  forces  become  nonlinear.  The 
governing  equations  were  reduced  to  ordinary  differential  equations 
by  Galerkin's  method  and  steady  state  solutions  were  determined  by 
harmonic  balance  technique. 

For  the  cantilever  beam,  the  solutions  agreed  qualitatively 
with  the  experimental  results  but  failed  to  predict  some  peculiar 
amplitude  drop-off  behavior  observed  in  the  experiment.  In  these 
cases  the  response  was  found  to  be  very  sensitive  to  small  changes 
in  initial  angle  of  attack. 

For  the  suspended  cable,  both  structural  nonlinearity  due 
to  tension  and  aerodynamic,  nonlinearity  had  to  be  considered.  The 
agreement  between  theory  and  experiment  was  good  for  the  frequency 
response  but  only  qualitative  for  the  amplitude  response. 
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SECTION  1 
INTRODUCTION 

Among  the  study  of  wind  effects  on  non-streamlined  struc¬ 
tures  (Refs.  1,  2,  3,  *0,  the  self-excited  oscillation  of  square 
section,  cylinder-like  bodies  has  been  a  subject  of  research  for 
many  years.  The  two  basic  exciting  forces  which  act  on  the  body 
are  due  to 

(a)  The  Harman- Vortex  shed  alternately  from  either 
side  which  exert  periodic  forces  on  the  body 

(b)  The  aerodynamic  lift  and  drag  forces  due  to 
asymmetric  flow,  including  the  effect  of  wake, 
flow  separation,  reattachment,  etc. 

The  vortex  shedding  is  responsible  for  high  frequency,  low 
amplitude  oscillations  often  called  Aeolian  vibrations  and  has 
been  observed  in  bluff  bodies  with  or  without  sharp  edges  such  as 
stretched  cables,  tall  towers  such  as  rockets  in  launching  pad, 
smoke  stacks,  etc.  The  oscillations  usually  occur  when  the 
natural  frequency  of  the  flexible  body  is  close  to  the  vortex 
shedding  frequency  (Ref.  5)*  Since  the  Strouhal  number  is  usually 
known  (Ref.  6,  7,  8),  this  kind  of  resonance  vibration  takes  place 
at  a  known  wind  velocity  and  can  easily  be  avoided. 

The  self-excitation  due  to  the  aerodynamic  lift  and  drag 
forces  usually  occurs  in  bluff  bodies  with  sharp  edges  and  is 
commonly  known  as  galloping.  It  has  been  observed  in  a  space 
shuttle  on  a  launching  pad  (Ref.  9),  tall  buildings  (Ref.  2)  and 
overhead  power  transmission  lines  during  a  sleet  storm  (Refs.  10- 
13).  This  kind  of  self-excitation  can  take  place  over  a  wide  range 
of  wind  velocity  independent  of  the  natural  frequency  of  the  body 
and  is  the  main  subject  of  study  in  the  present  report,  using  a 
square  section  vertical  cantilever  beam  and  a  suspended  sagging 
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cable  in  a  uniform  wind  stream  free  from  turbulence. 

The  aerodynamic  properties  of  bluff  bodies  exhibiting 
galloping  oscillation  are  such  that  at  a  certain  initial  angle  of 
attack  the  aerodynamic  lift  force  induced  by  a  small  transverse 
motion  of  the  body,  acts  in  the  direction  of  the  motion.  This 
results  in  a  negative  damping  type  exciting  force.  The  explana¬ 
tion  was  first  given  by  Den  Hartog  (Refs.  14,  15).  Parkinson  and 
associates  (Refs.  16-21)  have  analyzed  the  one  degree  of  freedom 
transverse  oscillation  of  square  section  prism  suspended  by  linear 
spring,  using  quasi-static  representation  of  aerodynamic  forces. 
The  resulting  nonlinear  equation  is  solved  by  Krylov  Bogoliubov’s 
asymptotic  method.  A  good  theory  to  experiment  agreement  is 
observed  except  when  the  structural  damping  is  very  small.  Novak 
applied  the  quasi-static  theory  using  energy  balance  to  a  two- 
dimensional  body  in  form  of  a  rigid  square  section  prism  pivoted 
at  the  base  (Refs.  22,  23)  and  obtained  qualitative  agreement  with 
the  experiment  and  did  not  get  the  hysteresis  in  the  experiment 
originally  observed  by  Parkinson.  Working  with  rectangular  pris¬ 
matic  pivoted  cantilever,  Novak  again  obtained  qualitative  agree¬ 
ment  with  the  experiment  but  observed  some  interesting  response 
characteristics  (Refs.  24,  25).  In  all  the  above  studies  the 
specimen  was  constrained  to  move  only  in  transverse  direction  and 
were  initially  kept  at  zero  angle  of  attack. 

In  the  present  report,  a  square  section  cantilever  beam  is 
considered  which  is  free  to  move  in  any  direction  without  twist¬ 
ing.  Quasi-static  representation  of  aerodynamic  forces  is  used 
to  determine  the  response  at  various  initial  angles  of  attack, 
theoretically.  Experiment  is  conducted  to  examine  the  applicabil¬ 
ity  of  such  a  theory  for  large  3D  displacement.  The  system  is 
governed  by  a  pair  of  nonlinear  partial  differential  equations. 

The  nonlinearity  arises  from  nonlinear  aerodynamic  forces.  They 
are  reduced  to  ordinary  differential  equations  by  Galerkin’s 
method  and  solved  by  harmonic  balance  technique  to  determine 
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steady-state  limit  cycle  amplitude.  Qualitative  agreement  with 
experiment  is  observed  at  low  amplitudes,  but  as  the  wind  veloc¬ 
ity  is  increased  peculiar  response  is  observed  in  certain  cases 
which  cannot  be  explained  by  quasi-static  theory.  The  static 
forces  are  found  to  be  quite  sensitive  to  slight  imperfections, 
sharpness  of  the  edges  and  Reynold's  number.  It  is  believed  that 
the  forces  are  even  more  sensitive  to  these  factors  when  the  body 
is  in  motion  and  the  averaged  static  force  measurements  are  not 
sufficient  to  explain  them.  The  quasi-static  representation  is 
explained  in  Section  2.  Sections  3  and  contain  the  cantilever 
analysis  and  experimental  results. 

The  galloping  of  power  transmission  cables  have  been 
studied  by  Hogg  (Refs.  11-12),  Richards  (Ref.  13),  and  many  other 
authors.  Richardson  and  Martuccelli  (Refs.  26-29)  worked  out  two- 
dimensional,  three  degrees  of  freedom  dynamics  of  a  D  section  with 
linear  spring  and  piecewise  linear  quasi-static  aerodynamics. 

They  also  measured  aerodynamic  properties  of  two-dimensional  sec¬ 
tions  of  various  shapes.  None  of  these  studies  consider  the 
aerodynamic  excitation  of  the  power  transmission  line  as  a  three- 
dimensional  catenary  in  writing  the  elastic  part  of  the  equations 
In  this  case  the  tension  in  the  cable  along  with  the  non-linear 
aerodynamic  forces  are  responsible  for  limiting  the  amplitude. 

Thus  structural  nonlinearity  arising  from  the  tension  terms  is 
also  important  in  this  case. 

In  the  present  report,  galloping  of  a  sagging  cable  of 
square  cross  section  has  been  obtained  experimentally  in  a  wind 
tunnel.  The  motion  is  governed  by  two  non-linear  partial  differ¬ 
ential  equations  which  are  reduced  to  non-linear  ordinary 
differential  equations  in  time  using  Galerkin's  approximation 
and  are  solved  by  harmonic  balance  method  using  only  two  equations. 
Good  agreement  is  obtained  in  the  frequency  of  oscillations  but 
the  amplitudes  are  in  qualitative  agreement  only  for  two  particu¬ 
lar  sags.  From  the  experimental  results  it  is  observed  that  the 


aerodynamic  lift  forces  do  not  reach  the  values  given  by  the 
quasi-static  theory  due  to  the  lack  of  spanwise  correlation.  The 
cable  analysis  and  experiment  are  described  in  Sections  5  and  6. 
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SECTION  2 


AERODYNAMIC  FORCES  ON  SQUARE  SECTION  BODY 


2.1  Introduction 

In  this  section  the  various  aspects  of  flow  aroung  a  fixed 
or  movable  square  section  body  and  the  nature  of  the  resulting 
forces  has  been  discussed.  Attention  is  restricted  to  uniform 
subsonic  flow.  The  flow  around  bluff  bodies  with  sharp  edges 
have  some  different  characteristics  in  comparison  to  the  flow 
around  a  cylinders.  These  are  pointed  out  in  Section  2,2.  The 
quasi-static  representation  of  aerodynamic  forces  on  a  translating 
square  section  body  is  explained  in  Section  2.3  and  is  used  in 
subsequent  sections. 

2.2  Aerodynamic  Forces 

Many  experimental  investigations  of  flow  around  a  circular 
section  body  have  been  done  .in  the  past  but  available  experimental 
results  for  a  square  section  body  are  relatively  few  even  though 
the  flow  around  bluff  bodies  with  sharp  edges  have  some  Interest¬ 
ing  characteristics.  The  basic  differences  between  the  two  flow 
fields  are  the  following; 

(a)  In  a  circular  body  the  flow  separation  can  occur 
at  different  locations  of  the  body  depending  on 
the  Reynolds  number  and  the  flow  condition  of  the 
wake.  In  square  body,  separation  can  occur  only 
at  sharp  edges . 

(b)  Because  of  fixed  separation  point,  there  is  a 
stronger  -spanwise  correlation  of  forces  on  a 
square  cylinder  than  on  a  round  one  (Ref.  33) 

Thus,  on  square  section  cylinders,  the  overall 
aerodynamic  forces  are  stronger. 


(c)  When  the  flow  is  at  an  angle  of  attack  due  to 
the  asymmetry  of  shear  layer  separation  from  the 
edges  and  presence  of  the  afterbody  in  the  sepa¬ 
rated  flow,  the  lift  forces  on  the  square  sec¬ 
tion  are  considerably  greater  than  that  for  the 
circular  cross  section  (Ref.  33). 

(d)  Reattachment  of  separated  flow  is  possible  at  a 
certain  angle  of  attack  and  separation  point  can 
shift  from  front  edge  1  to  the  rear  edge  2  at 
higher  angles  of  attack  as  shown  here  in  Pig.  3 
(Ref.  19).  Due  to  considerable  change  in  wake 
width,  the  drag  varies  considerably  with  angle 
of  attack  reaching  minimum  value  at  the  angle 

of  attack  when  reattachment  first  occurs 
(Pig.  3c  ).  The  lift  forces  reach  a  maxima  at 
this  angle. 

The  effect  of  the  Reynolds  number  on  a  square  section 
is  similar  to  that  of  circular  section  except  for  the  drop  in  drag 
at  critical  Reynolds  number  *1  x  10^  since  the  separation  point 
is  fixed.  However,  when  the  edges  are  slightly  rounded,  similar 
drop  in  drag  at  Re  6  x  10^  (for  r/b  =  0.176)  has  been  observed 
by  Delany  (Ref.  3*0 .  The  Strouhal  number  is  0.13  at  Re  =  3000 
and  changes  very  slightly  with  angle  of  attack  and  Reynolds  number 
(Refs.  6,  35). 

When  the  bluff  body  is  flexible  and  can  oscillate  trans¬ 
versely,  the  spanwise  correlation  of  forces  increases  even  with  a 
circular  cylinder.  The  shedded  vortex  can  excite  the  body  if  the 
shedding  frequency  is  near  its  natural  frequency.  The  oscillat¬ 
ing  body  in  turn  exerts  some  control  on  the  shedding  phenomena 
and  the  shedding  frequency  locks  on  to  the  oscillation  frequency 
over  a  short  range  of  velocity  instead  of  following  the  Strouhal 
number  relationship  (Refs.  5,  20,  36).  Extensive  study  of  wake 
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V, 


behind  an  oscillating  circular  cylinder  has  been  done  recently  by 
Toebes  (Ref.  37)  and  Griffin  (Ref.  38).  The  vortex  shedding 
frequency  of  the  present  specimen  is  quite  high  compared  to  its 
natural  frequency  and  thus  does  not  affect  the  response  appreci¬ 
ably. 

2.3  Quasi-Static  Representation 

Let  us  consider  a  square  section  body  translating  in  a 
uniform  fluid  stream  with  a  velocity  VJ  along  the  Z  direction  as 
shown  in  Fig.  2c.  The  induced  angle  of  attack  3  due  to  the  veloc¬ 
ity  is  given  by  tan (5  =  W/U.  In  quasi-static  representation  of 
aerodynamic  forces  it  is  assumed  that  the  forces  are  the  same  as 
if  the  body  were  at  rest'  in  a  uniform  stream  of  velocity  UREL  at 
an  angle  of  attack  3.  Total  forces  normal  to  and  along  the 

UREL  are 


L  -  57°UL  b  CL(f ) 


RE.L 

,2 


Ibs/fct 


(2.1) 


In  actual  flow  the  development  of  forces  due  to  an  induced  angle 
of  attack  may  not  be  instantaneous.  It  may  differ  both  in  magni¬ 
tude  and  phase.  In  quasi-static  representation  these  effects  are 
neglected.  Since  no  theory  exists  to  predict  the  unsteady  aero¬ 
dynamic  forces  on  a  bluff-body  with  turbulent  wake,  the  unsteady 
forces  can  only  be  determined  by  carefully  conducted  dynamic  tests 
For  square  section  body,  no  such  results  seem  available. 

Resolving  the  forces  along  Z  and  Y  direction  which  are 
normal  to  and  along  the  free  stream  direction  and  noting  that 
UREL  =  U  sees 

^2  "  g  f3  U2  b  Sec  p  |"  -  C\_  Cor>^i  —  Cx>  J 

Fy  =•  ^/°l)2b  -Seep  [  CDCosp  ~CLSmp]  1,\?) 
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Fg  is  taken  positive  in  the  direction  W.  To  use  Fz,  Fy  in  equa¬ 
tions  of  motion  it  is  convenient  to  express  them  as  algebraic 

• 

functions  of  W/U,  i.e.,  as  algebraic  functions  of  tan0.  One 
possible  way  is  to  define  new  force  coefficients  acting  in  the 
direction  Fz  and  Fy  such  that 

- 

j/u2b  CN(p) 

-c" 

11 

U2  k  Ccfp) 

(2.3) 

where  C^,  Cc  are 

defined  as 

cM(p)- 

Secji  [-Q(f)  - 

CEtan^  ] 

Cc(p)= 

Sec  P  [  cv(f)  - 

CLCcmp] 

(2.4) 

Considering  the  zero  reference  angle  of  attack  to  be  the 
position  when  one  of  the  faces  of  the  square  is  normal  to  the 
free  stream  direction,  C^,  Cc  with  respect  to  this  origin  can 
now  be  determined  as  following.  Lift  and  drag  forces  on  a  static 
body  at  various  angles  of  attack  a  (Fig.  2b)  can  be  measured  and 
converted  to  CL,  CD  by  dividing  by  ^/>U2b.  CN,  Cc  are  then 
computed  using  Eq.  2.4  .  Here  0  is  taken  equal  to  the  static 
angle  of  attack  a  at  which  the  forces  were  measured.  CN,  Cc  are 
now  plotted  against  tanji>  and  fitted  with  a  polynomial  in  the  form 

2  ( tan  f>  )n 

Y)=o  ' 

Z.  Qn(tanp)  (2.5) 

ms© 

In  the  equations  of  motion  C^,  are  used  as  force  coefficients 
and  tan0  is  replaced  by  W/U. 
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^r^T^T,  - 


•vaKycwpiwvc^vw^^  -~ 


Now  let  the  body  be  initially  at  an  angle  of  attack  aQ 

with  respect  to  the  free  stream  and  translate  without  rotation 

•  • 

with  velocity  components  V  and  W  along  Y  and  Z  directions  as 
shown  in  Fig.  2d.  Expressing  F2,  FY  in  terms  of  induced  angle 
of  attack  3  in  the  similar  manner  one  obtains 

Fi  =  b  [  -Cu(F-«.)Coof.  -  CD  (p-«<)  Sw.f>] 

Fy  =i/>U«Lb[  C„(p-T.)C=>,>|l  -  cL(p-«i)  5u»p]  (2'6) 


where 


a0  =  Fixed  angle  (+Ve  in  the  sense  shown  in  Fig.  2d) 
tan3  =  W/(U-V)  (Induced  angle  of  attack) 

UREL  =  ^  SecP 


(’)  =  > 


Here  CL(G-a0)  represents  static  CL  measured  at  an  angle  of  attack 
(G-a0)  where  zero  angle  of  attack  is  the  position  when  one  of  the 
faces  of  the  square  cross  section  is  normal  to  the  free  stream 
direction.  Introducing  UREL 

Fz  =f/3(U-V)Zb  [-Cl(P-Oo)  “^(p-cOto/np  ]Sec(b 

Fy  =^/°(U-V)2b  [  CD(p-&o)  -  CL(ji-oC)  tamp  ]Secf>  (2.7) 

»,  • 

Now  to  express  the  forces  in  terms  of  W  Find  V  in  a  convenient 
form,  it  has  to  be  written  in  a  polynomial  form  in  tanG.  This  can 
be  done  in  two  ways: 

a)  One  can  plot  [-CL(3-a0)  ~  C^(G-a0)tanG]SecG,  and 
[CD(3~a0)  -  CL(3-a0)tan3]SecG,  versus  tanG  at  a  fixed  a0  and  fit 
it  with  polynomials  of  tanG  to  give: 

_  N 


i/>u2b(i-g)  y  z>.)  tjp 


n=o 
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(2.8) 


Fy  =  ^/>uzb('-§f  I  Y.<°0  tcw>"p 

The  highest  order  power  N  can  be  chosen  for  a  best  reasonable 
fit  but  must  be  kept  minimum  to  avoid  unnecessary  complexity,, 

This  procedure  is  not  very  convenient  since  at  every  a0}  one  h-,.s 
to  carry  out  the  curve  fitting  numerically. 

b)  To  avoid  this  disadvantage,  it  is  desirable  to  express 
Eq.2.7  in  terms  of  and  Cr  whose  polynominal  fit  is  already 
known  from  Eq.  2.5.  To  do  this  one  may  write 

tom  p  -  (  tom  ( jb  -  cCi)  +  "to/vo  f  (  t  ~  tom  (  p  -  °£<s)  tcw>  0t,) 

(  />«.(/*-«•)  |-  tom(f>-tfo)  tom**)  (2.9) 

in  Eq.  2.7  and  combine  the  and  Cc  terms  using  the  definition 
(2.4)  of  CN  and  Cc  to  get 

2 

Fz  =  [Cn(P-M  -Cc(p-^to^o^]>>£cad/(l~tcw(p-cYe)  tow*.)2* 

FY  =  if  U  l  (»  -  §  f  [  C  c  ( f '  «.)  +  CM  ( f  -  «.)  U  - *«/(  I  -  MP- *)  tow  Ct.f 

(2.10) 


This  may  be  written  in  the  final  form  as 

Fz  0 -  %  f  [C»  (P“  «•)  ~  Cc.  (  p  ~  *o)  tom  Ofe  ]  ( l  +  tom  f  tomCC,f  Coo\a 

F  y  [pc  (p'tfo) -(-Cn(P“*o)  tom0[o^(  1+  tomptom«e,)a(^/33cfC) 

(2.xl) 
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Now  introducing  Eq.  2.5  in  above,  one  obtains 

f2  =  if <«.*«.  f  ?„ 

y^O  C  I  +  tcM\p>  t OMQo) 

Fy  -4 /> u*b ( I ■ - i  f CoJ«.  f  S„  1M: *?*«•>* 

u  ^  (l  +  Rw.pFoMCY.)"-1 


n?o 


(2.12) 


where  Pr>  =  PA  “  Qv\  bwcy0 

Q«  =  Qyy  +  Pnt'tVY'Ofo 


Eq.  2.12  can  be  expressed  in  the  form  (2.8)  as  follows.  First, 
one  may  write 


( 1  +  tcw\ 


I -0-2)  tbw|»to/»tf0  +  (tcwf>  to.-na'.y 


(2.13) 


For  example,  for  a  square  section  the  estimated  maximum 
I  tar\p  fcvnqfol  ^  0.05  .  Thus  If  N  =  6,  just  one  term  in  the  above 
approximation  needs  to  be  retained  with  only  3%  error.  Now 
using  Eq.  2.13  in  Eq.  2.12,  expanding  the  numerator  terms  fully 
and  collecting  the  coefficients  of  tan  jl  one  obtains 

F?  =  */>U,b(l-$)*£  Z.  {«.)  to-Tp 

n** 

F^^('4f?Y,K)U>  (2.„ 

w=c 

Here  Z  ,  Yn  are  known  functions  of  OC0 ,  Pn,  Qn  and  can  be  computed 
for  any  (X0  •  Z  ((Xo)>  Y  (<X0)  for  M  =  4  are  given  in  Appendix  C. 

XI  11  •  • 

In  the  equations  of  motion  tan^is  replaced  by  W/(U-V) 
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SECTION  3 


CANTILEVER  BEAM 

3.1  Introduction 

In  this  section  the  quasi-static  representation  of  the 
aerodynamic  forces  described  in  Section  2  will  be  used  to 
formulate  the  motion  of  a  square-cross  section  cantilever 
beam  when  placed  in  an  uniform  wind  stream.  Such  a  simple 
structure  is  chosen  so  that  the  emphasis  can  be  given  to  the 
quasi-static  representation  of  the  aerodynamic  forces  and  its 
applicability  can  be  examined  thoroughly.  The  geometric  non¬ 
linearity  due  to  large  deflection  of  a  cantilever  beam  is  very 
small  and  the  force-deflection  relation  is  linear  up  to  tip 
deflections  of  25%  of  the  span.  Thus  only  the  aerodynamic 
nonlinearity  needs  to  be  considered  here. 

The  beam  has  the  same  bending  rigidity  in  all  directions 
and  is  thus  free  to  translate  in  a  two  degrees  of  freedom 
space.  It  has  high  torsional  stiffness  and  does  not  twist 
under  the  existing  forces. 

Even  thoug.i  the  present  model  is  free  to  move  in  any 
direction,  the  motion  is  predominantly  normal  to  the  wind 
stream  direction  when  the  square  section  beam  is  initially 
at  zero  degree  angle  of  attack.  In  Section  3.2  this  one 
degree  of  freedom  motion  is  considered  and  solved  by  harmonic 
balance  method  for  initial  introduction  to  the  problem.  The 
solutions  are  compared  with  experimental  results  in  Section  4. 
This  is  the  only  case  investigated  by  most  authors.  Moreover, 
their  models  are  constrained  to  move  in  lateral  direction 
only.  In  Section  3.3  the  general  case  of  two  degrees  of  freedom 
motion  is  considered  and  solved  by  harmonic  balance  method.  The 
essential  features  are  pointed  out  with  a  simple  example  and  the 
justification  of  using  one  degree  of  freedom  motion  at  zero 
initial  angle  of  attack  is  discussed.  The  solutions  are 


1? 


compared  with  experimental  results  in  Section  4. 

3.2  One  Degree  of  Freedom  Motion 

Let  us  first  consider  the  case  when  the  beam  is  initially 
at  zero  angle  of  attack  with  one  of  the  sides  of  the  square 
section  facing  the  wind  stream  and  moves  normal  to  the  wind 
direction  as  shown  in  Fig.  2c.  The  simple-,  equation  of  motion 
of  the  uniform  beam  is  given  by  the.  well  known  expression 

mr&  +  cfx,  *  EI&  -  */•«**> o.(&/u)  (3-D 

Here  w(x,t)  is  the  displacement  normal  to  the  wind  stream,  U  is 
the  free  stream  velocity  and  c  is  small  structural  damping 
coefficient.  It  is  convenient  to  express  this  in  the  form  of 
viscous  damping  since  the  vibration  frequency  is  close  to  the 
natural  frequency.  In  general,  at  low  subsonic  speed,  the 
aerodynamic  forces  are  small  compared  to  the  inertia  and 
elastic  forces,  m  is  the  mass  per  unit  length,  p  is  the  den¬ 
sity  of  air  which  is  considered  incompressible,  b  is  the  dimen¬ 
sion  of  the  side  of  the  square,  and  the  right-hand  side  repre¬ 
sents  aerodynamic  .'orces  induced  by  the  motion.  Applying 
Galerkin's  approximation  to  reduce  Eq.  3.1  into  an  ordinary 
differential  equation  in  time  t  using  only  one  mode 


=  bfooyt) 


5  -  K/i 


(3.2) 


where  f ^  ( -; )  is  the  first  bending  mode  shape  of  a  freely  vibra¬ 
ting  cantilever  beam  as  given  in  Appendix  A,  the  Eq.  3.1 
becomes  after  nondimensionalization 


{.urn, 


(3.3) 
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where 


U)i  t  =  C 

-£^  (  First  mode  natural  frequency  ) 

A l  t t\(4 

4*,  =.  C  / 2(0^  (  Damping  ratio  ) 

^  =  y°ba/W  (  Mass  ratio  ) 

U.  =  U  f  6J,  b  (  Reduced  velocity  ) 


and  A^  and  are  integrals  of  the  mode  shape  and  are  defined 
in  Appendix  A. 

In  quasi-static  representation  of  aerodynamic  forces 

the  static  normal  force  coefficient  C„  versus  tan0  curve  is 

N 

represented  by  a  polynomial  in  tan8  in  the  form  Eq.  2.5 

cN=  i  p,(upf  =  <3-4> 

n-'O  T):0 

Introducing  Eq.  3,4  into  Eq.  3,3  nives 


%  +  X  +  X 


N 

^  Gn(%/u) 


Y\  =  0 


(3.5) 


1H 


where 


The  numerical  values  of  An  are  given  in  Appendix  A.  Subtracting 
out  the  static  part  from  above  and  using  to  denote  the  dynamic 
part,  Eq.  3.5  becomes 

N 

\  +  %\  -l/>u  2  G-(  Vu)  +  v°  (3-6) 

n-  z 

where 

%  =  (if,  -  jfv-P.) 

The  terms  under  summation  sign  contain  all  the  nonlinear  terms. 
The  solution  of  the  linear  equation  is  given  by 

Z  - - 

^  =  a  e  S‘/w-  (  J1~(<|>)  ^  +  f  )  (3.7) 

a  and  <|>  are  arbitrary  constants.  When  g^>0,  the  oscillations 
die  out.  Beyond  the  critical  velocity  uc=4s^/pP^,  when  g^  <0, 
the  oscillations  grow  and  full  nonlinear  equation  need  to  be 
considered.  If  the  amplitude  is  plotted  against  time  in  a 
semi  log  paper  from  an  experimental  transient  record,  the 
initial  slope  should  be  linear  and  equal  to  the  (^-puP^-s^)  if 
the  quasi-static  theory  works.  In  general  |g^|<<l  and  the 
frequency  of  oscillation  is  nearly  equal  to  the  natural  fre¬ 
quency. 

The  solution  of  full  Eq.  3.6  is  discussed  in  detail  by 
Bogoliubov  (Ref.  32) ,  Parkinson  (Ref.  18)  and  many  other  authors 
using  the  Krylov  and  Bogoliubov* s  method  of  slowly  varying 
amplitude  and  phase.  The  equations  can  have  stable  and  unstable 
limit  cycle  solutions,  depending  on  the  values  of  Gn.  They 
appear  as  concentric  circles  in  the  phase  plane  whose  origin 
is  an  unstable  focus.  In  this  report  Eq.  3.6  is  solved  by 
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harmonic  balance  method  to  determine  the  limit  cycle  solutions. 
The  stability  is  checked  by  slightly  perturbing  the  3teady 
state  solution.  All  these  methods  involve  the  assumption  that 
the  nonlinearity  is  small  compared  to  the  inertia  and  elasticity 
terms  and  the  oscillations  are  close  to  the  sinusoidal  motion 
near  the  natural  frequency.  For  a  square  section  body  in  low 
subsonic  flow,  these  are  realistic  assumptions. 

Let  the  steady  state  solution  be  given  by 


=  Q,  Si*  C  -V  b,  Co-6  ^ 


(3.8) 


Substituting  Eq.  3.8  in  Eq.  3.6,  collecting  the  coefficients  of 

sinx  and  cost  and  setting  them  to  zero  one  obtains 

Sint: 


- % b'  + 


xi— \ 
2.  .  ,  2  \  2. 


***,5,» 


b,(a?+b,) 


Cost: 


2  .  2 


<ja,  £2  a.(of+  b,1) 


*=3,5/ 


where 


g  =  G  .  C3)f0 

[  2  2-1 

Setting  a=  y  ,  Eq.  3.9  reduces  to, 


^  a  _  iy 'LL2  ^  Q,  (O'/u) 


n. 


=  0 


*0,5/ 


(3.9) 


(3.10) 


This  nonlinear  algebraic  equation  in  a  can  now  be  solved  to  give 
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the  limit  cycle  amplitude  at  a  given  u  and  4^.  Eq.  3.10  can 
be  expressed  in  another  nondimensional  form 


where 


I  s 

m  =  \ .  a  «  . 


=  o 


a  =  a 


U  =  u  C/r/4Ci) 


(3.11) 


For  a  given  (a/u) ,  one  can  calculate  u  without  solving  Eq.  3.11 
for  a.  The  resultant  plot  of  a  vs.  u  serves  as  a  master  plot 
for  a  particular  CN  vs.  tang  data  (3.4).  This  was  suggested 
by  Novak  (Ref.  23).  It  may  be  noted  that  (a/u)=(a/u)  and  the 
induced  angle  of  attack  6  at  the  tip  of  the  beam  is  given  by 
tan0=2 (a/u)  for  our  chosen  mode  shape.  Thus  the  steady  state 
solution  is  given  by 


%  =  f,(?)aS  mC 


(3.12) 


at  a  reduced  velocity  usu/w^b.  For  example,  for  a  fifth  degree 
polynomial  fit  in  Eq.  3.4,  Eq.  3.11  takes  the  form 

i(l)-[p'(l)  +  '•7S1  p>(tf  +  P5(|)5]  =  O 

Eq.  3.10  gives  both  stable  and  unstable  limit  cycle 
solutions  so  each  of  the  solutions  must  be  tested  for  stability 
Giving  a  small  perturbation  to  the  steady  state  solution,  ^  is 
written  in  the  form 

=  (Qs  +  S(t))  Sim  C  (3.13) 
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Here  a  is  a  steady  state  solution  satisfying  Eq.  3.10  and  £(t) 
s 

is  a  small  perturbation  from  the  steady  state  solution.  It  is 
also  assumed  to  be  a  slowly  varying  function  which  means  that 
the  increase  of  such  a  function  during  one  period  is  small  com¬ 
pared  to  the  average  value  for  that  period 


5  (C)  £  7T 

*Cc) 


«  I 


g  te)  2n 


«  1 


(3.14) 


According  to  this  assumption 


=  La»  +  5Cc)]  Co.  t 

oo 

“  ~[°i  +  ?CO]  Si«C  +  2£C z)C»,c 


(3.15) 


Substituting  this  in  Eq.  3.6  coefficients  of  sint  get  cancelled. 
Setting  coefficients  of  cost  to  zero,  cancelling  out  the  steady 
state  solution  and  keeping  only  the  linear  terms  in  £(t),  one 
obtains 

2|(r)  +  ®  ?(c)  =0  (3.16) 

CJ  Oj 

where 

F(o-)  =  2?,a  —  ^u2  £  (a/u)u 

5»- 

It  may  be  noted  that  F(a  )=0  gives  the  steady  state  solution. 

S 

Thus  in  order  that  the  steady  state  solution  be  stable,  one 
must  have  positive  at  F(a)=0. 

A  graphical  representation  of  the  steady  state  solution 
and  stability  is  explained  below.  Using  the  nondimensional 
notation  of  (3.11) 
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vf.  -  *seaBes.V'*^& 


F(a) 

dF(a) 

da 


-*/“>[*©  -  E  8-(f)  ] 


"  =  ‘,3,5,- 


l  i  -  &S5 !  I  =-©')  ] 

v>  M  -  I  4  e 


(3.17) 


Plotting  ^GnCa/u)  versus  (a/u)  ,  steady  state  solutions 
at  a  given  u  are  obtained  from  the  intersection  with  a  straight 
line  from  origin  and  with  slope  1/u.  The  solution  is  stable 
if  the  slope  of  the  straight  line  is  greater  than  the  slope  of 
the  curve  at  the  point  of  intersection. 

Gw(5/u)  plot  is  as  shown  in  Fig.  16a,  B2 
represents  unstable  limit  cycle  solution.  Fig.  16b  shows  the 
nature  of  steady  state  amplitude  variation  with  u.  The  dotted 
portion  represents  unstable  limit  cycle  solution.  At  large 
velocity  the  amplitude  is  asymptotic  to  a  line  of  slope  D. 

If  the  '^GnCo./v.)*  plot  is  as  shown  in  Fig.  16c  all  the 
steady  state  solutions  are  stable.  If  this  is  as  shown  in 
Figs.  16e  and  16g,  A2  represents  the  unstable  solution.  These 
different  forms  of  amplitude  response  for  bodies  of  various 
cross  section  have  been  discussed  by  Novak  (Ref.  24)  and  sta¬ 
bility  proven  by  the  method  of  Krylov  and  Bogoliubov.  By  the 
present  method,  the  steady  state  amplitude  and  stability  can 
be  determined  just  by  looking  at  the  (o-fv-)*  plot* 

3.3  Two  Degrees  of  Freedom  Motion 


ns  i,  v 


In  Section  3.2,  the  motion  of  the  vertical  cantilever 
beam  at  zero  initial  angle  of  attack  is  assumed  to  be  in  a 
plane  normal  to  the  free  stream  direction.  This  is  found  to 
be  a  reasonable  assumption  experimentally  for  this  cantilever 
beam  which  has  no  directional  preference  of  motion  structurally. 


in 


However,  when  the  beam  is  initially  at  an  angle  of  attack,  (see 
Fig.  2d) ,  the  motion  is  generally  in  a  plane  slightly  inclined 
to  the  z  direction.  This  inclination  is  presumably  decided  by 
the  aerodynamic  forces  alone.  Assuming  that  the  beam  can 
vibrate  freely  in  all  directions  with  the  same  natural  fre¬ 
quency  and  does  not  undergo  any  torsional  motion  under  the 
existing  forces,  the  equation  of  motion  can  be  written  in  the 
same  form  as  Sq.  3.1,  namely. 


n*  at  ax-i 


=  WbC2(f*/u  $/„  («„) 


(3.18) 


yr\  +  c  — 

a  t*  at 


4-  ELI 


=  WbcY(H/u,  Ufa,*.) 


Here  w(x,t)  is  the  displacement  normal  to  the  wind  direction; 
v(x,t)  is  the  displacement  along  the  wind  direction;  c  is  the 
structural  damping  coefficient  expressed  in  the  form  of  vis¬ 
cous  damping;  a0  is  the  initial  angle  of  attack  and  the  right 
hand  terms  represent  aerodynamic  forces  induced  by  the  motion. 
Again,  applying  Galerkin's  approximation  to  reduce  Eq.  3.18 
to  ordinary  differential  equation  in  time  t  using  only  one 
mode  for  w  and  v 


w(?,t)  =  bl.U^U) 
vf?.t)  =  bfumt) 


(3.19) 


5  =  X/j*. 


where  f^(£)  is  the  first  bending  mode  shape  as  in  Eq.  3.2,  the 
Eq.  3.18  becomes  after  some  nondimensionalization 

i‘t2f  ?  +  q, 


%’ 2  s  i  m,  -  ( I#>i  ■  V-)  l «« 5 


(3.20) 
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where  GJ,  t  =  C 

=  c/ 2Yn<a, 

(°)  s=  dO/dC 

Other  notations  are  the  same  as  in  Eq.  3.3.  In  Section  2.3  a 

quasi-static  representation  of  the  motion  induced  aerodynamic 

•  • 

forces  on  a  body  translating  with  velocity  W  and  V  without 

rotation  has  been  derived.  Using  th*  Eq.  2.14,  C  and  Cv  can 

2  Y 

be  expressed  in  the  form 

cH  =  0~£)2!  (M>' 


>2  £ 


(3.21) 


CY  s  ^ 


for  small  tan  Stance. 


($f(- 


.  .-(cm -z) 

V) 

vj 


(3.22) 


In  general,  because  of  presence  of  damping  due  to  drag,  the 
stream-wise  motion  is  very  small,  and  one  may  write 

o-*r- 1 + (^(^) + -  (3.23) 

« 

In  the  present  problem  |V/u|<0.1.  Thus  if  M=7,  only  two  terms 
need  to  be  retained  in  the  above  approximation  with  only  4% 
error  . C z,  C can  now  be  expressed  in  the  fc.m 


C2  =  2°  (J  ~  u  )  +  2‘ +  Z2-(§F 


(3.24) 


7  Z 


•  /  • .  n. 


1 1  ^  uum 


mi  =  3  =  o 


?1 


Cy  expression  is  the  same  as  C2  except  Zm  is  replaced  by  Ym* 

In  this  case  Nmo=l,  Nmf=(m-2)  N  =(m-2) (m-l)/2.  Eq.  3.24 
^  •  , 

expresses  C2  #  Cy  in  terms  of^,^  in  a  convenient  form  which 
can  be  converted  to  ordinary  different  equation  in  time  using 
Galerkin's  approximation.  Introducing  Eq.  3.24  in  Eq.  3.20/ 
one  obtains  after  subtracting  out  the  static  part 


VU+  %A-i? 


Vl 


mi*)  n 


Y  win  = 

Aw+m  — 
%  = 


Yrr,  Ni  wn  A  m+n 


W+n-H 


x,]Lf,(?)3  m 

O 

(2?,  -  ^uZ,) 


^12  =  (P'U-  H.) 

%>  =  (4/uY) 

til  =  (2?,  +  /^Yo) 


2? 


Here  q^  q2  denote  only  dynamic  terms  since  static  parts  have 
been  subtracted  out.  The  terms  in  the  square  bracket  contain 
all  the  nonlinear  terms.  The  mass  and  stiffness  matrix  are 
symmetric  and  positive  definite.  Thus  the  instability  can 
arise  only  from  the  damping  terms.  It  should  be  remembered 
that  the  self-excitation  is  mainly  due  to  negative  g^  and  not 
because  of  interaction  between  q^  and  q^. 

First  consider  the  linear  part  of  the  Eq.  3.25. 


+  ^  =  o 

&|  dll  b|  0|2.  o| 

a  +  +  +  9/  =o 

*2.  on  b|  dti  h  o  2. 


Let  the  solution  be  j=a^e^T ,  c||2=a2e^T*  Substituting  them  in 
Eq.  3.«26  a  pair  of  homogeneous  simultaneous  equations  in  a^ 


and  a2  are  obtained. 


0^+  %)>  +  1 ) 


?12>  a. 

o  K 


=  o 


(3.27) 


For  non- trivial  solution,  the  determinant  of  the  coefficients 
must  be  zero.  The  variable  parameter  is  the  reduced  velocity  u. 
It  may  be  noted  that  for  a  perfectly  symmetric  cross  section 
body  at  ao=0  g12=0  ^21=®  anc^  linear  equations  are  uncoupled. 
Since  g22  is  positive,  the  motion  in  streamwise  direction  die 
down.  Thus  if  g^  is  negative,  the  motion  normal  to  the  free 
stream  direction  is  only  excited.  The  equations  become  coupled 
if  the  square  section  is  ?+■  an  initial  angle  of  attack  aQ,  or 
if  P  and  Q.  are  not  zero  .  a  =0  (for  example,  if  the  section 
is  not  perfectly  symmetric) ,  From  experiment,  it  is  found 
that  the  motion  for  the  initial  angle  cases  is  not  normal  to 
the  wind  direction  but  in  a  plane  inclined  to  it.  Setting  the 
determinant  of  coefficients  of  Eq.  3*27  to  zero  a  fourth 


degree  polynomial  in  p  is  obtained  which  can  be  solved  to  give 


=r  <r  ±  i  J I  -  cr2-’ 


where 


(3.28) 


It  can  be  seen  that  at  low  velocity  when  (9;q922~^12-21^ 
both  the  values  of  a  are  negative  since  in  general  (g-^+g^) 
is  positive  for  a  square  cross  section  body,  at  all  velocity. 

The  oscillations  die  down.  As  u  is  increased,  flutter  condition 
is  reached  when 


(%  %i%> )  0 


(3.29) 


which  makes  cr=0  and  p=+i.  Thus  the  flutter  frequency  is  the 
same  as  the  natural  frequency  w-p  It  should,  however,  be 
remembered  that  the  self-excitation  occurs  because  of  the 
negative  damping  term  g^  and  not  due  to  interaction  of 
and  . 

The  Eq.  3.29  can  be  solved  for  u  to  give  the  critical 
velocity  u 

(2Y.-Z,)  +  /(2Y.-Z,)1  +  8(ZrY.  -Y,Z.) 


(2.Y.-Y,  H.) 


(3.30) 


ignoring  the  negative  value. 

At  u>uc  one  of  the  values  of  o  becomes  positive  and  the 
oscillations  grow.  The  motion  corresponding  to  the  negative 
o  die  out.  In  general,  | o | < < 1  and  the  frequency  of  oscillation 
is  close  to  the  natural  frequency.  Thus 


1° 


24 


0"  ±  1 


(3.31) 


Substituting  this  in  Eq. 


3.27  and  using  the  assumption  a 


2 


<<1 


0_2 

a, 


^21 

2crJ 


±  i8 

e 


t awn  6  ®-<r  (3.32) 


For  the  specimen  used  in  the  experiment,  at  low  subsonic  speed 
o  is  of  order  0.01  so  the  phase  difference  between  q^  and  q2 
at  low  amplitude  is  small  and  the  initial  motion  is  in  an 
inclined  plane.  The  phase  difference  between  q^  and  q2  repre¬ 
sents  an  elliptic  whirling  motion.  It  should  be  noticed  that 
at  critical  velocity  u^q^  and  q2  are  in  phase  and  the  motion 
is  purely  planar. 

Let  us  now  look  for  steady  state  solution  of  the  total 
Eq.  3.25.  Since  the  elasticity  terms  are  linear,  and  the 
linear  and  nonlinear  damping  terms  are  small,  the  frequency 
of  oscillation  beyond  the  critical  velocity  are  close  to  the 
natural  frequency  of  oscillation.  Thus  a  steady  state  solu¬ 
tion  to  equation  (3.25)  can  be  taken  as 

Q;  =  Q.i  iiw'  C  +■  Con  Z. 

11  (3.3 

^2.  =  Stvyi  C  +  b2£aoC 

Substituting  these  into  Eq.  3.25,  collecting  coefficients  of 
sinx ,  cost  and  setting  them  to  zero  one  obtains  four  algebraic 
equations  coming  only  from  the  first  derivative  terms. 

Si at  z  : 

Con  C  : 

Sin  Z  : 

Cot j  Z  : 

?r> 


-  j/V-f51(u,,b,,o2j  b2)  =  0 

~3icbl  -?ub2  =  O 

V'  *fnaz  =o 


(3.34  a,b,c,d) 


Here  f  ^ (a^/b^,a2 ,b2)  etc.  contain  all  the  nonlinear  terms  in 
al/bi/Ci/^i  an^  are  ?iven  Appendix  B  for  a  general  poly¬ 
nomial  fit  of  aerodynamic  data  for  the  case  b^=0.  In  steady 
state  limit  cycle  solution/  the  starting  point  of  time  is 
irrelevant  and  can  be  set  to  zero.  Setting  also  b2=0 
Eg.  3.34a/C  are  automatically  satisfied  and  one  is  left 
with  two  equations  (3,34b,d)  with  two  unknowns  a^  and  a2- 


%Q|  +  V*  =° 


(3.35  a,b) 


which  can  be  solved  at  a  given  u  to  obtain  the  steady  state 
solution. 


S  =  f,(5)  a,  Z»,z 


V 

b 


{,  (0  a2  Swn  £ 


(3.36) 


The  inclination  of  the  plane  of  the  motion  with  respect  tc  the 
z  axis  is  given  by 


=  tana.  (  0,2./  a. 


(3.37) 


Positive  y  has  the  same  sign  convention  as  positive  aQ.  The 
total  tip  amplitude  W.j,  is  given  by 


^  =  2at/Co^y 

b 


(3.38) 


It  may  be  noted  that  Eq.  3.35  can  be  nondimensionalized  in  the 
special  form  of  Eq.  3.11. 

To  examine  the  stability  of  the  steady  state  solution, 
we  again  assume 


2( 


(3.39) 


^  =  L  ais  +  He)]  Si«c 
\  -  [  a2.s  +  T  col  Si«c 


where  a^g,  a2s  are  steady  state  solutions  satisfying  Eq.  3.35 
and  C(t),  n  (t)  are  slowly  varying  functions  as  defined  in  Eq. 
3.14.  Thus  proceeding  in  the  same  manner 


X 

o 

X 


t  tt|»  +  SCO  ]  ConC 

[°M  +  ’He)  ]  GwC 


oo 

\  ~-[Qis  +  5(c) ]  sv*c  +2$tz)Co»e 

^  ^-[a25+  YlC^)]  SLnC  2  »l(c)  Cg/j  C 

(3.40) 

Substituting  these  in  Eq.  3.35,  the  coefficients  of  sinr  get 
cancelled.  Setting  coefficients  of  cost  to  zero,  cancelling 
out  the  steady  state  solution  and  keeping  only  the  linear  terms 
in  £(t)  and  ii(t),  one  obtains 


where 


l  O 
o  I 

a,t  = 

®  12.  — 


$ 

o 


a„  a 

0-2.1  Q 


11  J 


a 


21  = 


a2Z  = 


z  ( %  -  ll: ) 

a  t/i 


i  c%*~  MlO 

•^az 


=  0 


AT 


(3.41) 


a,  =  als 

Oz  —  C»2S 
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This  gives  the  familiar  eigenvalue  problem.  For  stability , 
the  root  of  the  characteristic  equation 


2. 

|°  +  \>  (  atl  +  ^22.)  +  (  atl°22-'  av2Q2l)  =  0  (3.42) 

must  have  -ve  real  part.  The  condition  is  givun  by 

(a„  Qn)  >  o 

(Q»Qn-Qi»a»<)  >  o  |3-43) 


Let  us  now 
are  given  by 


consider  a  simple  example  for  which  CN  and  CQ 
P(  Tam  |5  +  to/n^  |i> 

Q0  +  tom  ji  -V-  (3.44, 


In  order  to  get  a  feel  for  the  quantities  involved,  P^O.9, 
P3=-8,  Qq=1.24,  Q2=-3,  Q^=30  roughly  fits  the  experimental  data 
shown  in  Figs.  7a,  8a  up  to  tans=0.45  if  one  ignores  the  small 
off-symmetries.  In  one  degree  of  freedom  case  this  form  of 
CN  has  the  response  characteristic  as  shown  in  Fig.  16c  in 
Section  3.2.  There  is  one  stable  limit  cycle  solution  and  the 
limit  cycle  amplitude  increases  monotonically  with  u.  With  the 
limitation  j tan6tanaQ | <0 . 1 ,  |V/U | <0 . 1  the  harmonic  balance  equa¬ 
tions  (3.35)  are 


5„a' +  Hzf +  8.m  (&) 

+  4., 7  Zs(i)5  ] 

+  Vl~2?*£l7C'Y>(u)3+  8  34 

+  Z3&)5] 


-  O 


(3.45) 


=  0 
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with  the  approximate  quasi-static  representation.  Adopting 
the  nondimensional  form  of  (3.11),  Eq.  3.45  become  after 
division  by  (a,/Ui) 


(i-z.)  +  2Z4  - 


rn  .2 


+  1.17  Z5(|f  J  =  0 


(3.46  a,b) 


-v,  +  ( i + 


r-['-7«'Y»(|)a  +  8.31 

+ 4-,7Y5(irj  =  o 

=  Xt(  //4£) 

=  a,C//4<,) 

=  uy  =  f 


’  s  \  4 


where 


OX.  = 


Qi  _ 

a, 


Here  Z  .  .Zc ,  \  .  .Yc  are  functions  of  P  .  P.  ,  P.,,  Q  .  Q~,  Q. 
and  aQ  and  are  given  in  Appendix  C.  Since  CN  is  cubic,  the 
second  and  third  terms  in  square  bracket  of  Eq.  3.46a  are 
very  small  compared  to  the  first  term  and  can  be  neglected. 
From  Eq.  3.46b 

x/  .  .  v/  (M1  .  .  v/  /a,\4 


~  Y,+  4  1.t7Y5(|) 

torn  y  -  y  -  /-  s  4 

l+2*)-8.3lY,(£) 


(3.47) 


In  determining  (a^/u)  by  substituting  (3.47)  in  (3.46a),  one 
may  ignore  the  fourth  power  terms  and  get 


I  +  2^/(1 +  2Y,) 

J  '-7«'  {  Z3-  2  Z„Y5/(  i  +  2Y„)  ] 


The  limit  cycle  solution  is  possible  only  if  the  term  under  the 
radical  sign  is  positive.  If  CN  is  antisymmetric  and  Cc 
symmetric  in  6  then  at  aQ=0  y^,  y^,  y5  are  zero  and  the  solu¬ 
tion  is  given  by 


£  =  M7^ 0/wT 

fowl  y  =  o 


(3.49) 


If  04=0,  one  may  ignore  y3,  y$.  Then  tany  is  approximately 

given  by 

tew  y  ~  7—^ —  ~  ~  U-  |i)fei^0Ce  (3.50) 

—  +  2Y  z  To  v  Uo  / 

U  0 

Thus  qualitatively  if  no  drag  variation  with  $  is  considered 
then  y saQ;  if  the  drag  increases  with  S,  Q2>0,  Y<ao;  the 

drag  decreases  with  8,  Q2<0,  '*’>a0*  When  Q4  is  included  the 
problem  is  more  involved.  Numerical  examples  are  considered 
in  Section  4.  In  general,  y  does  not  affect  a^  very  much,  a^ 
is  primarily  determined  by  Z^  and  E^. 
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SECTION  4 


EXPERIMENT  WITH  CANTILEVER  BEAN 


4.1  Introduction 

Since  no  purely  theoretical  method  of  predicting  the  aero 
dynamic  forces  on  a  blunt  body  with  a  turbulent  wake  exists, 
experimental  observations  form  the  basis  of  the  present  investi¬ 
gation.  The  quasi-static  representation  of  the  aerodynamic 
forces  constitute  one  of  the  means  of  predicting  the  response  of 
a  self-excited  square  section  body.  So  experiments  have  ben 
conducted  with  a  square  section  cantilever  beam  placed  in  a  uni¬ 
form  wind  stream.  First  the  static  lift  and  drag  forces  are 
determined  and  are  used  in  quasi-static  theory  to  determine  the 
responses.  Then  tests  are  carried  out  to  determine  the  response 
experimentally  with  the  specimen  at  various  initial  angles  of 
attack.  Some  interesting  response  phenomena  are  observed.  They 
are  described  in  Section  4.3.  In  Section  4.4  the  results  obtained 
by  this  quasi-static  theory  using  the  static  test  lift  and  drag 
data  are  compared  with  experiment.  Only  qualitative  agreement 
has  been  observed. 

A  movie  of  the  experiment  has  been  taken,  which  shows  the 
self-excitation  phenomena.  The  transient  motion  and  the  response 
at  various  angles  of  attack  can  be  clearly  observed. 

4.2  Static  Force  Measurement 

The  experiment  with  the  cantilever  beam  is  carried  out  in 
a  1  ft.  x  1  ft.  low  speed  wind  tunnel.  The  air  is  blown  into  the 
tunnel  by  a  blower  and  is  smoothed  out  by  a  series  of  netting  and 
honeycomb  mesh.  The  velocity  is  increased  by  the  area  contraction 
ratio  of  36:1.  The  flow  is  smooth  and  nearly  free  from  turbulence. 
The  specimen  is  placed  in  the  open  jet  1.5  inches  from  the  end  of 
tunnel  as  shown  in  Figs.  1  and  2a.  All  the  velocities  mentioned 
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are  the  velocities  just  inside  the  tunnel  and  are  used  as  the 
reference  velocity  throughout  the  report.  Before  the  testing 
the  velocity  near  the  specimen  position  is  measured  with  a  pitot 
tube  whose  total  pressure  holes  are  at  1.5  inches  from  the  end 
of  the  tunnel  and  static  pressure  holes  1"  farther  behind.  In 
Fig.  JJ,  this  is  plotted  against  the  reference  velocity.  The  re¬ 
duction  in  velocity  at  the  specimen  position  due  to  the  open  jet 
effect  is  found  to  be  less  than  2%  at  all  reference  speeds  up  to 
80  ft/sec.  and  can  be  considered  within  tolerable  experimental 
error. 

The  specimen  cantilever  beam  is  made  of  1/8  inch  diameter 
steel  rod  and  is  13.1  inches  long.  It  is  coated  with  wax  and 
drawn  through  a  square  template  to  give  it  a  square  cross  section 
outside.  The  side  of  the  square  is  0.09  inches.  The  section  is 
not  perfectly  square  but  is  slightly  rounded  at  the  edges.  The 
wax  coating  is  thin  and  does  not  change  the  structural  property 
of  the  rod  appreciably.  Thus  the  rod  can  be  assumed  to  have  the 
same  area  moment  of  inertia  in  all  directions.  It  is  free  to 
deflect  in  any  direction  without  twisting  and  with  some  bending 
stiffness.  This  simplifies  the  theoretical  analysis  considerably. 
The  13.1  inches  long  specimen  is  clamped  in  a  vertical  position 
and  the  base  is  bolted  to  the  top  of  a  turn  table  which  can  be 
rotated  through  360°.  The  orientation  of  the  specimen  with  respect 
to  the  wind  stream  direction  and  the  chosen  coordinate  axes  are 
shown  in  Figs.  2a,  2b.  The  four  sides  are  identified  as  side 
A,  B,  C,  D  with  respect  to  a  small  120  ohm  strain  gage  which  is 
placed  at  the  root  of  the  beam  on  side  B  to  record  the  transient 
motion.  When  the  side  A  is  facing  the  wind  stream,  the  general 
position  is  termed  as  position  A  and  so  on.  The  static  and  dy¬ 
namic  calibration  of  the  strain  gage  is  shown  in  Fig.  9  at 
position  A,  aQ  =  0  for  different  tip  deflection  WT  in  Z  direction, 
and  is  quite  linear.  There  are  small  variations  if  the  motion  is 
slightly  inclined  from  Z  direction. 


The  clamped  end  of  the  specimen  is  2  inches  below  the  wind 
tunnel  floor  to  avoid  flow  disturbance  from  the  base  support  fix¬ 
tures.  Since  there  is  little  motion  at  the  root,  the  motion  in¬ 
duced  forces  are  negligible  and  the  set  up  does  not  affect  the 
response  significantly.  Tests  are  first  conducted  to  determine 
the  static  lift  and  drag  forces  at  various  angles  of  attack.  This 
is  done  by  measuring  the  static  tip  deflection  in  Z  and  Y  direc¬ 
tions  visually,  at  different  angle  of  attack  position  a  as  shown 
in  Fig.  2b,  and  then  calculating  the  force,  assuming  uniform  load 
over  the  span  of  the  beam  exposed  to  the  wind,  using  linear  force 
deflection  relation. 

CL  =  8E1  Mr/ if  I  -  1(f)1  +  >  (f  )«  ] 

C„  =  8ElV7/4y.u'^['  -|(f)5+i(fn 

(4.1) 

Here  A£  is  the  portion  of  the  beam  at  the  root  which  is 
outside  the  wind  stream.  For  AJl  =  2  inches,  J,  =  13.1  inches, 
the  factor  [  I  -  0.9955,  showing  that  the 

effect  is  very  small.  C^,  CD  are  computed  at  37,  44,  56,  65  ft/sec 
velocity  for  position  A  and  at  37,  45  ft/sec  velocity  for  position 
D,  and  are  shown  in  Figs.  5  and  6.  The  variation  in  C^,  CD  from 
the  error  in  measurement  of  deflection  is  less  than  5$.  It  is 
observed  that  while  remains  relatively  unaffected  by  velocity, 
there  appears  a  Reynolds  number  effect  on  CD.  For  an  increase  in 
Reynolds  number  from  1800  to  3100,  CD  increases  from  1.24  to  1.54 
at  a  =  0.  The  square  section  CL  is  characterized  by  the  negative 
lift  curve  slope  in  the  vicinity  of  zero  degree  angle  of  attack. 

The  lift  attains  peak  value  near  o  =  12°  where  drag  reaches  a 
minimum.  In  this  vicinity  the  flow  first  reattaches  on  one  side 
and  the  wake  has  minimum  width.  Beyond  this  angle  the  flow  is 
fully  attached  to  both  the  windward  faces  but  the  wake  width 
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increases  as  the  square  is  rotated  farther,  and  the  drag  increases 
simultaneously.  Even  though  the  general  nature  of  C^,  variation 
with  a  are  similar  for  the  position  A  and  D,  the  differences  are 
noticeable  at  small  angles  of  attack,  signifying  that  the  static 
forces  are  quite  sensitive  to  the  imperfections  and  sharpness  of 
the  edges. 

Using  the  C^,  CD  data,  (^(g),  Cc(g)  as  defined  in  Eq.  2.4 
are  computed  and  plotted  against  tannin  Figs.  7  and  8.  An  8th 
degree  polynomial  least  square  fit  is  shown  in  solid  line.  The 
fit  is  good  only  up  to  tang  =  +  0.45  and  does  not  fit  properly 
near  sharp  peaks.  The  coefficients  of  the  polynomials  Pn  and  Qn 
are  given  in  Table  2. 

4.3  Description  of  Experimental  Result 

The  first  natural  bending  frequency  of  the  beam  is  9.4 
cycles/sec.  The  theoretical  estimation  with  bending  rigidity  for 
a  1/8  inch  round  steel  rod  gives  9.5  cy/sec  showing  that  the  wax 
coating  does  not  affect  the  structural  stiffness.  First  the 
critical  damping  ratio  is  obtained  from  the  dynamic  decay  record 
by  plotting  the  transient  amplitude  against  time  on  semi-log 
paper,  and  is  shown  in  Fig.  10  along  with  a  typical  decay  record. 
The  damping  ratio  is  constant  at  very  low  amplitude  but  varies 
linearly  at  higher  amplitude.  The  constant  value  is  due  to  vis¬ 
cous  and  structural  damping.  The  linear  variation  is  due  to  the 
aerodynamic  damping  as  the  body  moves  through  still  air.  If  it 
is  assumed  that  the  drag  at  any  point  of  the  harmonically  oscillat¬ 
ing  beam  is  proportional  to  the  square  of  the  velocity  at  that 
point,  then  it  can  be  shown  by  equating  the  work  done  over  a  cycle, 
the  equivalent  damping  ratio  is  proportional  to  the  amplitude  of 
the  cycle. 
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(4.2) 


WT  is  the  tip  amplitude  and  p",  A^,  A2,  w,  are  as  defined  in 
Eq.  3.3.  This  is  plotted  in  Pig.  10b  from  the  origin  for  =  1.6 
and  2.  Since  e<<l,  o;/u>,  *»  1  ,  it  can  be  seen  that  the  slopes  of 
linear  variation  are  much  smaller  than  the  experimental  plot. 
For  CD  =  1.6,  the  SEa  =»  O.87  x  10 “\wT/t>) whereas  the  slope  of  ex¬ 
perimental  linear  variation  £ea  »  3.8  x  10”\w-r/b)which  is  four 
times  as  much.  It  is  not  quite  clear  why  this  is  so  but  one  of 
the  reasons  could  be  that  the  beam  is  actually  not  moving  through 
still  air  but  sets  the  wind  in  motion  as  it  moves  through  it.  As 
it  swings  back,  it  meets  the  incoming  turbulent  wind  which  in 
effect  doubles  the  net  velocity  and  the  drag  is  four  times  as 
much.  It  should  be  noted  that  when  the  beam  is  oscillating  lat¬ 
erally  in  a  wind  stream  this  air  damping  in  that  direction  is  no 
longer  present.  Only  the  structural  damping  which  is  obtained 
from  the  intercept  of  the  experimental  curve  on  the  <5,  axis,  impedes 
the  motion  in  this  situation.  This  is  estimated  to  be  between 
.002  and  .0025  in  terms  of  critical  damping  ratio.  The  value  .002 
is  used  for  all  calculations. 

The  response  at  zero  degree  initial  angle  of  attack  is 
measured  for  all  the  four  positions  A,  B,  C,  D,  i.e.,  with  each 
of  the  four  sides  facing  the  wind  stream.  The  turn  table  is 
clamped  at  the  position  A  at  zero  degree  angle  of  attack  and  tunnel 
is  turned  on.  At  a  speed  below  the  critical  velocity,  there  is 
no  response.  When  the  beam  is  disturbed,  the  oscillations  die 
down.  The  wind  velocity  is  then  increased  in  small  steps.  At 
around  10  ft/sec,  the  beam  begins  to  oscillate  by  itself.  The 
amplitude  increases  slowly  and  then  settles  down  to  a  steady  state. 
To  record  the  transient  motion  at  a  particular  velocity,  the  beam 
is  brought  to  rest  by  touching  it  at  the  top  and  then  letting  go. 
The  Sanborn  recorder  which  is  hooked  up  to  the  strain  gage  is 
turned  on  simultaneously.  The  amplitude  grows  slowly  and  then 
settles  down  to  a  steady  state  limit  cycle.  The  final  amplitude 
is  quite  steady,  signifying  absence  of  large  scale  turbulence  in 
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the  flow.  The  beam  is  then  given  a  small  kick  to  see  if  it  comes 
back  to  the  same  steady  state  value.  It  is  then  given  a  big  kick 
to  examine  if  there  are  other  limit  cycle  solutions  at  that  veloc¬ 
ity.  The  steady  state  tip  amplitude  is  measured  visually  with  a 
scale  held  on  top  of  the  tunnel  outside  the  wind  stream.  Some  of 
these  measurements  are  compared  with  static  calibration  of  the 
strain  gage  done  at  the  beginning  of  the  experiment  for  various 
tip  deflections  due  to  an  end  load,  and  is  given  in  Fig.  9.  They 
are  almost  the  same  if  motion  is  nearly  normal  to  the  wind  stream 
direction.  The  steady  state  amplitude  and  transient  response  are 
recorded  for  the  position  A  for  velocities  up  to  65  ft/sec.  The 
velocity  is  then  reduced  in  steps  and  the  amplitudes  are  measured. 
The  process  is  repeated  at  zero  degree  angle  of  attack  with  sides 

B,  C,  D  facing  the  wind  stream,  but  no  transient  records  are  taken. 
The  result  is  plotted  in  Fig.  11a.  WT  is  the  tip  amplitude.  It 

is  observed  that  even  for  amplitudes  as  large  as  WT/b  =  31,  the 
frequency  is  still  nearly  9.4  cycles/second  as  in  small  amplitude 
oscillation.  The  steady  state  amplitude  variation  with  velocity 
is  more  or  less  the  same  for  all  the  four  sides  up  to  about 
25  ft/sec.  Beyond  this  velocity  there  is  a  marked  departure.  In 
position  A,  the  amplitude  levels  off  and  then  drops  as  the  veloc¬ 
ity  is  increased  to  45  ft/sec  but  rises  again  as  the  velocity  is 
increased  farther.  Similar  response  variation  is  observed  for 
the  position  B  and  C  but  the  drop  off  and  rise  occurs  at  differ¬ 
ent  velocities.  For  position  D,  however,  the  amplitude  rises 
monotonically  with  velocity  to  a  large  amplitude  of  WT/b  =  31  at 
55  ft/sec.  The  specimen  is  not  tested  beyond  this  velocity  to 
avoid  its  destruction.  With  two  dimensional  square  section  body, 
response  behavior  of  form  shown  in  Figs.  16a,  16c  have  been 
observed  but  this  type  of  behavior  is  rather  unusual.  In  order 
to  investigate  this  peculiar  amplitude  response  in  position  A,  B, 

C,  further  tests  are  conducted  for  the  specimen  position  A  and  D 
so  that  the  basic  differences  can  be  detected. 
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Now  instead  of  keeping  the  specimen  at  zero  initial  angle  of 
of  attack,  it  is  placed  at  a  small  initial  angle  of  attack c t0  as 
shown  in  Fib .  3b  for  each  of  the  positions  A  and  D  and  the  steady 
state  amplitudes  are  measured  at  various  velocities.  Now,  however, 
the  motion  is  no  longer  in  the  Z  direction  but  generally  in  a 
plane  slightly  inclined  to  it.  The  total  tip  amplitude  W,„  and 
the  inclination  of  the  plane  of  the  motion  with  z  direction  denoted 
by  V  ,  are  plotted  in  Figs.  12a-f  (position  A)  and  12g,h  (position 
D)  for  various  initial  angles  of  attack  position  oc0  and  at  differ¬ 
ent  velocities.  The  sign  convention  for  i  and  ctD  are  as  shown 
in  inset  diagram  of  Fig.  12a.  They  are  both  positive  in  the  same 
sense.  At  position  A,  for  velocities  below  25  ft/sec,  the  steady 
state  amplitudes  are  nearly  symmetric  for  positive  and  negative 
0io  and  goes  down  as  0£o  is  increased  (Fig.  12a).  It  goes  down  to 
zero  ir.  the  vicinity  of  =  10°.  AtOC0  beyond  this  the  oscilla¬ 
tions  die  out.  The  plane  of  oscillation  is  however  somewhat 
irregular.  At  17  ft/sec,  the  /  is  in  opposite  sense  to(X0.  At 
24  ft/sec  the  inclination  \  are  nearly  zero.  However  beyond  this 
velocity,  there  is  a  marked  departure  from  this  pattern  as  shown 
in  Figs.  12c-f.  The  responses  are  no  longer  symmetric  with  re¬ 
spect  to  o cQ  and  there  is  a  large  change  in  amplitude  and  inclina¬ 
tion  of  plane  of  motion  for  a  small  change  in  from  -2.5°  to 
+2.5°.  The  peak  amplitudes  occur  in  the  vicinity  of  (X,Q  =  +7.5° 
and  goes  through  a  minima  near  aQ  =  -2.5°  (Figs.  12c, d).  There 
is  also  considerable  change  in  inclination  of  plane  of  oscillation 
)(.  At  small  &0j  K  is  in  the  same  sense  as  0to  but  bigger.  Beyond 
WQ  =  +7.5°  as  the  amplitudes  drop  of f , y  has  the  tendency  of  going 
down  to  zero  again.  As  velocity  is  increased  the  response  becomes 
more  sensitive  to  W0 .  At  65  ft/s  (Figs.  12e,f)  the  amplitude  has 
peaks  at  w0  =  2.5°,  +10°  and  drops  off  in  between  near  o£0  =  +5°. 
Near  <X0  =  +10°,  the  amplitude  and  inclination  are  rather  large  but 
drops  off  quickly  as  oc0  increased  a  little  farther.  The  same 
tests  are  carried  out  with  the  specimen  in  position  D,  at  37  and  45 
ft/sec  (Figs.  12g,h).  The  amplitude  in  this  case  shows  the 
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regular  trend  as  in  Fig.  12a,  with  peak  at  a0  =0,  dropping  off 
as  ac  is  increased.  Oscillation  ceases  near  aQ  =  +7°.  The  in¬ 
clinations  are  nearly  zero  but  show  trends  similar  to  those  of 
Fig.  12a.  Thus  it  can  be  observed  that  the  drop  off  and  rise  of 
amplitude  at  ct0  =  0  for  position  A  (Fig.  11a)  is  associated  with 
large  sensitivity  of  the  response  for  small  change  in  a0  in  the 
vicinity  aQ  =  0,  whereas  in  the  portion  where  the  amplitude  rises 
monotonically  with  increase  in  U,  a  small  change  in  a0  from  a0  =  0 
does  not  change  the  amplitude  or  inclination  considrably.  The 
amplituue  drops  off  quickly  as  a0  is  increased  farther. 

Fig.  13  shews  some  sample  transient  records  for  position  A, 
a0  =  0.  They  are  plotted  in  sami-log  paper  on  real  time  scale  in 
Fig.  I1!.  The  initial  logarithmic  increment  in  amplitude  is  found 
to  be  linear.  The  experimental  values  of  initial  linear  log 
(Slope)/03.  at  different  free  stream  velocity  are  plotted  in  Fig. 

15  with  a  negative  sign.  According  to  quasi-static  theory  this 
is  the  experiment?1  value  of  g^/2,  the  linear  theory  -  equivalent 
-  damping  coeficient  as  defined  in  Eq.  3.6. 


All  the  experimental  results  described  here  have  been  found 
to  be  quite  repeatable. 


H.*l  Theory  ar.d  Experiment  Comparison 


Using  the  static  aerodynamic  force  measurements  it  is  now 
possible  to  calculate  the  response  characteristics  using  the 
method  described  in  Section  3>  and  compare  with  the  experiment. 
When  the  beam  is  initially  at  zero  angle  of  attack,  the  steady 
state  amplitude  at  a  given  velocity  is  given  by: 


i(i\ . 

u  lu  j 


p#) +  p>i)1+  /U7 psd)s+  »•«  fvd)71 


=  0 


(*1.3) 

for  a  cantilever  beam  with  the  mode  shape  used  and  with  a  8th 
degree  polynomial  fit  of  versus  tan3  curve.  Using  the  8th  de¬ 
gree  polynomial  fit  of  versus  tanB  of 
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a)  Fig.  7a,  tested  at  37  ft/sec  in  position  A 

b)  Fig.  7b,  tested  at  37  ft/sec  in  position  D 

the  steady  state  amplitudes  a  are  plotted  against  v.  in  Fig.  ?.lb, 

along  with  the  experimental  values  obtained  for  the  specimen 

in  position  A  and  D  at  ocQ  =  0.  The  value  £-j  =  .002  p>  -  0.35x10"^ 
has  been  used.  Pn  are  given  in  Table  2.  It  should  be  noted 
that  2(a/-u)  represent  the  induced  angle  of  attack  ji>  at  the  tip 
of  the  beam.  Since  the  cui've  fitting  is  only  up  to  tan  /*  =  0.45, 
results  only  up  toS/u=  0.225  are  valid.  The  results  plotted  are 
wit.iin  this  range.  For  a  quick  conversion  from  physical  units 
to  nondimensional  values,  the  following  relations  may  be  used  for 
0.04  38 

U  O-'T/SEC)  =  O.M3U  =  10.1  u 


Wj 

b 


2  a  =  ^5.6  a 


(4.4) 


Reasonably  good  agreement  is  observed  for  amplitudes  up  to 
WT/b  s  15  .  The  experimental  amplitude  variation  in  position  D 
shows  similar  trend  but  are  substantially  lower  in  magnitude.  The 
quasi-static  theory  fails  to  predict  the  drop-off  characteristics 
ooserved  in  position  A  and  in  general  over  estimates  the  steady 
state  amplitude. 

Another  way  of  examining  the  applicability  of  quasi-static 
theory  is  to  compare  the  initial  transient  growth  rate  when  the 
amplitudes  are  small  and  equations  are  linear.  To  do  this,  the 
linear  theory-equivalent-damping  ratio  p^/2  =  o.25y®u  P,  ,  as 

defined  in  Eq.  3*6  is  plotted  in  Fig.  15  and  compared  with 
experimental  values.  P^  is  the  initial  slope  of  Cf>J  curve  (Fig.  7a 
U  =  37F/S,  Fig.  7b)  and  =  .002,  They  are  of  the  same  order  of 
magnitude  and  trend  but  not  close  to  linear  with  velocity  as  given 
by  the  theory.  It  may  thus  be  concluded  that  the  actual  aero¬ 
dynamic  forces  are  more  complex  than  a  simple  quasi-steady  descrip 
tion.  The  static  lift  and  drag  values  do  not  contain  all  the 
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information  required  to  describe  the  dynamic  forces  induced  by  the 
motion.  It  is  probable  that  as  the  amplitude  increases  beyond 
WT/b  =  10,  the  flow  is  no  longer  two  dimensional  and  a  break  in 
the  flow  pattern  at  a  part  of  the  beam  causes  a  flow  change  in 
the  entire  portion. 

The  amplitude  response  at  different  a 0  can  be  compared  with 
the  theory  in  a  qualitative  basis  since  in  the  one  degree  of 
freedom  case  it  has  been  observed  that  the  quasi-steady  theory 
cannot  predict  accurately  the  amplitude  response  when  the  ampli¬ 
tudes  are  large,  but  gives  the  right  trend.  So  an  approximate 
fit  of  and  Cc  data  has  been  chosen  for  comparison  with  the 
experiment  (see  Pigs.  73  8). 

CN  =  0.  9  ta/Yi  jb  -  8.0taY\5|b 
Cc  =  1.24  -  3.0  +  30.0 

CN  =  1.2  tcwjb  -  I5.0fa/n3(3 
Cc  =  >•'  -  3.0  ta/v^p  +  30.0  p 

After  a  steady  state  solution  has  been  obtained  at  a  given  a 0  from 
Eqs.  3.47,  3.48  one  should  check  that  +  !  taA\oro|  is  less 

than  the  limits  of  curve  fit  (4.5).,  lta*p  tcw«r<>|MA<o.  |  t  (v/l7l  <  o.| 
in  order  that  the  solution  be  valid.  The  curve  fit  is  not  very 
good,  but  it  is  observed  that  the  subsequent  value  of  /  does  not 
affect  the  amplitudes  substantially.  The  value  of  V  is  however 
quite  sensitive  to  Cc«  The  theoretical  results  are  plotted  in 
Fig.  12  in  broken  lines  for  comparison  with  experimental  values. 
Results  at  velocity  17.9  ft/sec  is  riot  computed  since  the  Cc  data 
is  not  known.  Amplitudes  at  velocities  24.5  and  37  ft/sec, 
position  A  (Figs.  12b, c)  and  37  and  45  ft/sec,  position  D 


(PCS.  A  ,  (lev*  p|  <  0.45) 


(POS.  13  ,  l  tcvYi p  |  <  04  ) 
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(Fig.  12g,h)  agree  with  the  experiment  qualitatively  but  the  in¬ 
clinations  do  not.  It  seems  a  simple  quadratic  of  Eq.  3*50 
can  explain  qualitatively  the  experimental  behavior  at  small  0Co 
since  it  gives  values  proportional  to  0i!o.  At  low  velocities 
the  dynamic  C’c  possibly  increases  with  j?>  which  makes  jf  go  in 
opposite  sense  too<s0.  At  higher  velocity  the  dynamic  Cc  possibly 
goes  down  with  fb  following  the  initial  static  trend  (Fig.  8a) 
making  Q2  negative  and  ^  goes  up  linearly  with  <XQ  for  small  0lo 
(Fig.  12c).  Same  explanation  apply  to  Fig.  12h  but  not  to  12g. 
Thus,  the  aerodynamic  drag  forces  acting  on  the  oscillating  body 
seem  to  be  rather  complex  and  do  not  necessarily  follow  the 
static  curve. 

From  these  comparison  it  is  observed  that  in  two  degrees 
of  freedom,  motion  the  quasi-static  theory  can  give  a  qualitative 
description  of  the  steady  state  amplitude  response  if  it  gives 
the  right  results  in  one  degree  of  freedom  case.  The  peculiar 
response  obtained  in  position  A  (Fig.  12d,  e,  f)  cannot  be  ex¬ 
plained  by  introducing  the  two-dimensional  motion. 

From  the  experiment  and  quasi-static  theory  comparison  it 
may  be  concluded  that  the  self-excitation  phenomena  of  a  square 
section  beam  can  be  described  qualitatively  by  the  theory  but 
usually  results  in  overestimation.  Physically  the  beam  cannot 
achieve  the  high  amplitudes  predicted  by  theory  without  going 
into  geometric  nonlinearity  but  below  the  velocity  30  ft/sec 
(u  =  3j  Fig.  lib)  it  can.  Up  to  W^./b  =  15  (a  =  0.3)  the  flow  is 
fairly  two-dimensional  and  the  amplitudes  follow  the  quasi-static 
results  fairly  well.  As  the  amplitude  increases,  the  top  portion 
of  the  beam  moves  in  an  arc  which  could  be  responsible  for  differ¬ 
ent  flow  patterns  and  lower  amplitudes.  On  a  bluff  body  with 
sharp  edges  where  the  flow  is  separated  and  the  wake  is  turbulent 
the  aerodynamic  forces  are  apparently  quite  complex  and  sensitive 
to  imperfections  and  sharoness  of  the  edges,  and  can  give  rise  to 
unpredictable  responses. 
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SECTION  5 


SUSPENDED  CABLE 

5.1  Introduction 

The  galloping  phenomena  of  overhead  suspended  power 
cables  has  been  a  subject  of  research  for  many  years.  The 
self-excitation  is  believed  to  be  due  to  the  negative  damping 
characteristic  of  the  cable  cross  section.  However,  unlike 
a  cantilever  beam  the  amplitudes  can  be  limited  by  both 

(a)  Nonlinear  aerodynamic  force  characteristics  of 
the  cross  section. 

(b)  The  tension  developed  in  the  beam  due  tc  large 
deflections . 

To  study  this,  galloping  of  a  sagging  cable  of  square  cross 
section  has  been  obtained  experimentally  in  a  wind  tunnel 
and  is  described  in  Section  6.  Theoretical  means  of  describing 
the  motion  using  quasi-static  aerodynamic  theory  is  discussed 
here.  The  motion  is  governed  by  two  nonlinear  partial  dif¬ 
ferential  equations  which  are  reduced  to  nonlinear  ordinary 
differential  equations  in  time  using  Galerkin's  approximation 
and  are  solved  by  harmonic  balance  method  using  only  two 
equations.  The  theory-experiment  comparison  is  discussed  in 
Section  6. 

5.2  Equations  of  Motion 

Consider  an  uniform  cable  which  is  suspended  between  two 
poles  with  a  particular  sag  and  placed  in  a  wind  stream  (Figs. 
17a,  b) .  With  the  wind  off,  the  cable  is  in  a  catenary  shape 
and  the  two  ends  are  clamped  in  that  position.  The  flexible 
cable  is  of  square  cross  section  and  in  the  suspended  position 
one  of  the  flat  sides  face  the  wind  stream.  It  is  assumed 
that  the  cable  can  deflect  in  y  and  z  direction  without  twisting 
(Fig.  I7c)  like  the  giant  wheel  riding  boxes.  Thus  only  two 
degrees  of  freedom  need  to  be  considered  in  the  equilibrium 
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equations.  Only  small  sags  up  to  midspan  sag/span  ratio  d/Jl 
of  0.08  are  considered.  Thus  the  slopes  at  a  deflected  posi¬ 
tion  are  small  and  the  simple  beam  equations  are  adequate  to 
describe  the  cable  motion.  Also,  for  small  sag  it  can  be 
assumed  that  the  tension  in  the  cable  is  constant  along  its 
whole  length  and  is  thus  independent  of  coordinate  x.  It  can 
be  shown  using  catenary  equations  that  for  a  simply  supported 
cable  without  bending  rigidity  and  under  constant  self  weight 
and  drag  forces  along  the  span,  the  tension  T* at  the  support 
and  at  midspan  are  given  by 
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sd 


'  +  (W™£)2  [  |  +  4(d/e)2] 
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(5  1) 


Here  m  is  the  mass  per  unit  length,  Z  is  the  span,  d  is  the  sag 
at  midspan.  DQ  is  the  drag  force  per  unit  length.  For 
d/£<.080,  the  total  variation  of  tension  along  the  span  is 
less  than  5%.  For  a  cable  with  small  bending  rigidity  it  is 
believed  that  the  assumption  of  constant  tension  does  not 
involve  error  more  than  this. 

When  the  wind  is  turned  on,  at  low  velocity  the  cable 
leans  back  and  takes  up  a  stable  static  position  until  the 
critical  velocity  i.*>  reached.  As  the  wind  velocity  is  increased 
beyond  the  critical  velocity,  the  oscillations  grow  until  a 
stable  limit  cycle  is  reached.  Thus  to  analyze  the  motion 
it  is  convenient  to  subtract  out  the  static  part  from  the 
overall  equation 
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Here  F2  ,  FY  represent  the  aerodynamic  forces  induced  by  the 
motion.  All  the  static  terms  of  aerodynamic  forces  are  in¬ 
cluded  in  mg  and  Dq  for  static  analysis.  Subtracting  the 
static  equations 
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from  Eq.  5.2,  one  obtains 
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The  Wo(x),  Vq(x)  satisfy  the  static  equation  (5.3)  along  with 
the  boundary  conditions 

AT  X  ~  O  ]  Wo  =  0 


dWo 

=  dWo 

dX 

dX 

(5.5) 
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where  W  (x)  is  the  deflection  shape  with  the  wind  off,  i.e., 
oo 
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Dq=0.  The  double  zero  subscript  will  be  used  to  refer  to 
quantities  when  the  free  stream  velocity  U=0.  The  W  and  V 
satisfy  clamped-clamped  boundary  conditions  at  both  ends.  The 
bending  rigidity  El  is  assumed  to  be  the  same  for  both  W  and 
V  displacements.  The  coupling  between  W  and  V  arises  from  the 
aerodynamic  and  tension  terms  only.  It  is  assumed  that  even 
after  U  exceeds  the  critical  velocity,  the  static  equations 
(5.3)  are  satisfied  and  give  a  measure  of  the  mean  position, 
i.e.,  the  position  if  no  motion  induced  forces  were  present. 
Since  the  tension  is  assumed  to  be  independent  of  x,  at  a 
fixed  time  t,  the  tension  T^  can  be  obtained  by  calculating 
the  strain  in  the  form 


Tt  =  Too  +  —  j  (6  5  -  d/i) 
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This  tension  term  is  the  source  of  nonlinearity  in  the  cable 
vibration  equation  apart  from  the  aerodynamic  nonlinear  terms, 
Both  of  these  are  responsible  for  limiting  the  amplitudes. 
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In  quasi-static  representation  of  aerodynamic  forces 


it  will  be  assumed  that  CN  and  Cc  are  in  the  form 


where 


C  N  Pjiour'P  +  P^  tew? 
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It  is  possible  to  work  out  the  problem  with  higher  order  poly¬ 
nomials  for  better  representation  of  CN  and  Cc  but  it  would 
result  in  added  complexity  without  adding  better  insight  to 
the  problem.  Proceeding  in  the  same  manner  as  in  Section  3.3 
F2  and  Fy  can  be  written  as 


f2  -  +  p3(i)5(,  +  £)] 

Fy  =  [-2^  +  ($)*  ]Q. 


(5.9) 


for  | V /y  |  <0.1.  Actual  values  of  \J/\j  are  still  smaller  and  it 
is  possible  to  neglect  the  underlined  terms  within  the  present 
approximation. 

5.3  Solution  of  Static  Equations 

The  static  equations  (5.3)  along  with  the  boundary  con¬ 
dition  (5.5)  can  be  solved  exactly  if  Tq  were  known.  In  order 
to  determine  Tq  it  has  to  be  expressed  in  terms  of  the  dis¬ 
placements  in  the  form  (5.7).  This  makes  Eq.  5.3  nonlinear  and 
can  only  be  solved  approximately. 

Solving  Eq.  5.3  in  terms  of  T  with  the  boundary  con¬ 
dition  (5.5) ,  one  obtains 
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The  double  zero  subscript  refers  to  quantities  when  the  wind 

is  off.  For  these  cables,  e  <e  <0.1.  Thus 

o  oo 

Using  this  approximation  Eq.  5.10  simplifies  to 
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The  above  is  valid  for  0<5<.5,  and  the  functions  are  assumed 

symmetric  on  the  other  half  of  the  span.  Now  initial  tension 

T  can  be  determined  by  setting  W  =d  at  £=0.5  and  then  solv- 
oo  J  v  oo 

ing  for  Tqo  with  the  above  approximation,  to  get 
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The  effect  of  bending  rigidity  is  to  reduce  the  tension  from 

2 

the  average  'Catenary  analysis' value  of  mg£  /8d.  In  the  present 

3 

case  EI/mg£  =0.01,  so  the  correction  is  of  the  same  order  as 
the  approximation  involved  in  the  analysis.  Now  the  static 
tension  Tq  at  a  particular  velocity  U  can  be  determined  by 
substituting  Eqs.  5.11  in  the  expression  for  Tq  in  Eq.  5.7  and 


carrying  out  the  integration.  A  nonlinear  algebraic  equation 
in  T0/T00  is  obtained  which  can  be  solved  by  iteration.  Sub¬ 
stituting  them  back  in  Eq.  5.11  would  give  the  static  deflec¬ 
tion  shapes  Wq(£),  Vq(5)* 

The  variation  of  T  /T  with  wind  velocity  U  is  in 

o  oo 

general  as  shown  in  Fig.  19.  For  a  simply-supported  cable 
with  no  bending  rigidity  the  problem  is  statically  determinate 
and  can  be  solved  exactly  using  the  'Catenary  equations’  to 


give 
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(5.13) 


This  is  compared  with  the  present  results  in  Fig.  19  and  is 
found  to  give  higher  values.  Using  Eg.  5.11  one  can  obtain 
the  static  positions  at  different  velocity.  One  such  plot 
for  d/£=0.059,  0^=1. 4  is  shown  in  Figa  20.  For  a  simply- 
supported  cable  the  static  positions  are  given  by  straight 
lines  inclined  to  the  vertical  at  an  angle  tan-1 (Q/mg) .  The 
inclinations  are  found  to  be  close  to  the  present  case  but  the 
static  positions  are  not. 

5.4  Harmonic  balance  solution 

Using  Galerkin's  approximation  Eq.  5.4  is  reduced  to 
ordinary  differential  equations  in  time  assuming  W  and  V  in  the 


form 
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Here  f^(5),  f ^ (£)  are  first  and  second  bending  mode  shape  of 
a  freely  vibrating  clamp-clamped  beam  and  are  given  in  Appendix 
D.  fQ(C)/  fQ(5)  are  static  deflection  shapes  given  by  Eq.  5.11. 
It  may  be  noted  that  f^(£),  fQ(£)  are  symmetric  functions 

about  midspan  while  f 2 (£)  is  an  anti-symmetric  function.  Under 
the  existing  forces  the  four  modes  are  sufficient  to  describe 
the  motion.  Hore  (Ref.  31)  have  analyzed  the  free  vibration  of 
a  similar  flexible  sagging  cable  using  the  present  mode  shapes 
and  have  obtained  good  theory-experiment  agreement. 


In  nondimensional  form,  the  ordinary  differential  equa¬ 
tions  can  be  expressed  as 
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■  The  time  is  nondimensional  with  some  reference  frequency  (o„. 
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The  definition  of  p,  u  are  slightly  different  than  those  used 

f  in  Section  3.  A  ,  B. , ,  B__,  D, ,  D~ ,  8  are  constants  and  are 

mn  lx  22  1  2  r 

given  in  Appendix  D.  qQ,  qQ  are  a  measure  of  static  position 
of  the  cable.  They  vary  with  velocity  as  shown  in  Fig.  21. 

Their  definition  are  also  given  in  Appendix  D.  It  may  be 
noted  that  r^,  r2  also  vary  with  velocity. 

So  far  no  structural  damping  terms  have  been  included. 
There  is  a  certain  amount  of  uncertainty  in  its  actual  values. 
Experimentally  with  the  wind  off  the  cable  can  be  given  small 
tapping  in  first  or  second  horizontal  bending  mode  and  the 
amplitude  decay  recorded.  The  slope  of  the  transient  amplitude 
plotted  against  time  in  a  semi-log  paper  is  found  to  be  linear. 
From  this  viscous  damping  coefficients  c^,  c2  can  be  evaluated. 

j  Assuming  c^c^,  c2=c2  the  structural  damping  can  be  expressed 

as  a  viscous  damping  for  theoretical  convenience  by  adding 
the  terms  (c,/wcj*)c^  }  (c,  /w  >  (c2  /m  uR  )  ^  ate.  to 

Eqs.  5.15..  However,  with  the  wind  on,  the  cable  will  oscillate 
at  a  higher  frequency  and  it  is  not  known  how  much  is  the  actual 
structural  damping  then.  To  overcome  this  difficulty  c^,  c^ 
etc.  can  later  be  set  to  zero  and  the  resulting  change  in 
responses  compared  since  the  structural  damping  are  usually 
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small.  Let  us  denote 
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The  first  and  second  mode  equations  decouple.  The  stiffness 
matrix  of  (5.17  a,b)  is  symmetric  and  positive  definite.  So 
the  instability  can  only  arise  from  negative  damping  terms 
g^  and  g2* 

Considering  Eq.  (5.17  a,b) ,  the  critical  velocity  and 
the  frequency  of  oscillation  at  the  critical  velocity  can  be 
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determined  in  the  usual  manner  by  introducing  c^=o,  e.  3 
and  setting  the  condition  of  condition  of  existence  of  non¬ 
trivial  solution  of  the  resulting  homogeneous  algebraic  equa¬ 
tions.  This  gives 


uoc 


03, 


~Z„ 

-ft- 1  ^  I 


at  a  velocity  which  makes 


.2-7  rs^  2.  /s^ 

+  n 


ftf-lL?)1  nff,  -t 

i  *,  +  *. 


=  0 


(5.19) 


(5.20) 


Since /l,  /  JL,,  ^  are  functions  cf  velocity  (Fig.  21),  no 

explicit  expression  for  critical  velocity  can  be  obtained,  and 
ujp,  uc  has  to  be  determined  by  plotting  (5.20)  for  various  u 
until  it  becomes  zero.  Corresponding  Wp  will  give  the  flutter 
frequency. 

The  Eq.  (5.17  c,d)  are  uncoupled  from  each  other.  The 
critical  velocity  for  ^  mode  will  be  given  by  g2=0.  T^e 
flutter  frequency  is  r2*  Since  g2  is  positive  at  all  velocity 
will  die  out. 

For  the  solution  of  entire  non-linear  equation  let  as 
take  only  Eqs.  5.15  a,b  assuming  no  second  mode  is  present. 

The  equation  has  quadratic  and  cubic  nonlinear  terms  in  cj,  and 
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The  reference  frequency  is  chosen  equal  to  the  linear  first 
mode  natural  frequency  at  u=0,  i.e. 
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So  at  u=0  r^l  in  Eq.  5.15. 

To  get  an  idea  of  the  nature  of  the  solution  of  Eq.  5.21 

is  first  integrated  numerically  using  predictor  corrector 
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scheme  with  an  error  bound  10  .  Some  results  are  shown  rn 

Fig.  25.  For  small  initial  displacements  the  amplitudes  grow 
initially  then  settle  down  to  a  steady  state  amplitude.  An 
estimation  of  induced  angle  of  attack  shows  that  the  induced 
angle  of  attack  (3  is  well  within  the  linear  range  of  CN, 
according  to  the  quasi-static  theory.  Thus  to  evaluate  the 
steady  state  amplitude  by  harmonic  balance  it  is  sufficient 
to  assume  (Ref.  33) 
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w  is  the  nondimensional  frequency  with  respect  to  ujr.  Sub¬ 
stituting  (5.23)  in  Eq.  (5.21),  keeping  only  the  coswt ,  sinwr, 
cos2wt,  sin2wT  terms  and  setting  their  coefficients  and  the 
constant  terms  to  zero  one  obtains  10  nonlinear  algebraic  equa¬ 
tions  in  12  unknowns  bQ,  b^...^  and  w,u.  Since  in  steady 
state  solution  absolute  phase  is  unimportant,  a^  is  set  to  zero 
and  equations  are  solved  by  Newton-Raphson  iteration  procedure 
for  the  10  unknowns  w,  bQ,  b^.b^  fc>2 ,  a2  at  a  given  u.  The 
algebraic  equations  and  details  of  iteration  procedure  are 
given  in  Appendix  E.  Some  of  the  solutions  are  shown  in 
Fig.  24  and  are  discussed  in  Section  6.3.  The  solution  of 
the  form  (5.23)  represents  large  swinging  motion  of  the  cable. 
Since  the  cable  swings  way  up  above  the  X-Y  plane  (Fig.  17  b,c)  , 


cos2ux  terms  come  in  for  representing  the  motion  in  the  q^, 
cartesian  coordinate.  From  the  polar  coordinate  point  of 
view  the  frequency  of  oscillation  of  the  physical  cable  is  ok 

Assuming  the  quasi-static  representation  of  the  motion 
induced  lift  forces  are  correct,  theoretically  the  amplitude 
is  limited  by  the  tension  together  with  the  q^  damping  term 
which  becomes  important  when  the  cable  moves  against  the  wind 
stream  at  the  end  of  the  larg^  swinging  motion. 


SECTION  6 


EXPERIMENT  WITH  SUSPENDED  CABLE 


6.1  Introduction 

To  analyze  the  galloping  phenomena  of  suspended  power 
cable,  tests  are  first  conducted  with  a  uniform  sagging  cable  of 
square  cross-section  in  a  wind  tunnel.  The  cable  is  flexible 
but  with  a  definite  but  small  bending  rigidity.  Large  and 
rather  violent  self-excited  oscillations  are  observed.  When  a 
cable  of  circular  cross-section  is  placed  in  the  wind  tunnel  no 
such  self-excitation  is  obtained.  Thus  the  galloping  phenomena 
is  found  to  be  due  to  negative  damping  type  aerodynamic  forces 
typically  associated  with  bluff  bodies  with  sharp  edges.  The 
cable  begins  to  oscillate  by  itself  beyond  a  certain  critical 
velocity  but  finally  settles  down  to  a  steady  state  amplitude  . 

The  experiments  are  conducted  with  different  sag/span  ratios  of 
d/i  =  0.080,  0.059,  0.043,  0.029,  0.020  and  0.011  and  the  results 
are  described  in  Section  6.2.  Theoretical  solutions  using  simple 
quasi-static  theory  are  compared  with  the  experimental  results 
In  Section  6.3. 

6.2  Description  of  Experimental  Results 

The  experiment  is  conducted  in  a  5  ft  x  7  ft  closed  circuit 
low-speed  wind  tunnel.  The  specimen  is  made  from  a  0.045  inch 
diameter  flexible  cable.  It  is  coated  with  wax  and  drawn  through 
a  template  with  a  square  hole  to  give  it  a  0.09  inch  square  sec¬ 
tion  outside.  Care  is  taken  to  avoid  twisting  of  the  cable  during 
the  drawing  process.  Because  of  this  construction  procedure  the 
uniformity  is  not  perfect  and  the  edges  are  somewhat  rounded.  The 
extension  rigidity  EA  of  the  bare  cable  has  been  measured  by 
Ilore  (Ref.  31)  and  is  used  here  assuming  that  the  wax  coating  does 
not  affect  it  appreciably.  The  bending  rigidity  HI  is  determined  by 
measuring  the  tip  deflection  due  to  self-weight  of  a  portion  of  the 


cable  cantilevered  at  one  end  and  using  the  linear  force  deflec¬ 
tion  relation  El  =  vn ft4  /  g  g  where  6  is  the  measured  tip  de¬ 
flection.  There  is  certain  amount  of  scatter  and  the  mean  value 
is  taken.  The  numerical  values  of  the  cable  parameters  are  given 
in  Table  3. 

The  cable  is  suspended  in  the  wind  tunnel  from  two  rigid 
1.5  feet  long  0.75  inch  diameter  poles  fixed  to  the  floor  of  the 
tunnel  (Pig.  17a).  The  cable  is  suspended  at  a  particular  sag 
and  is  clamped  at  two  ends  in  such  a  manner  that  the  cable 
assumes  a  natural  catenary  shape  and  one  of  the  sides  of  the 
square  faces  the  wind  stream.  Because  of  flexibility  and  clamped 
edge  condition  the  cable  can  deflect  in  Y  and  Z  directions  (Pigs. 
17b, c)  without  any  apparent  rotation  like  the  riding  boxes  of  a 
giant  wheel.  The  posts  are  long  enough  so  that  the  entire  cable 
is  outside  the  wind  tunnel  boundary  layer.  There  is  some  inter¬ 
ference  at  the  ends  because  of  the  clamps  and  the  posts.  At 
the  beginning  of  a  test  run  at  a  particular  sag,  the  first  hori¬ 
zontal  mode  natural  frev.aency  of  the  cable  is  measured  by  slightly 
tapping  the  cable.  This  frequency  gives  the  experimental  value 
of  linear  first  mode  natural  frequency  at  U  =0  as  given  in 
Eq.  5.22  and  is  denoted  by 

At  several  different  sag  positions  decay  "Q<  ,.-rds  for  small 
horizontal  oscillations  are  obtained  using  a  small  strain  gage 
at  one  end  of  the  cable.  The  strain  ga^e  signals  for  small  mid¬ 
span  deflections  are  linear.  From  '  uc  lope  of  the  transient 
amplitude  plotted  against  time  on  *.  •;  -log  paper,  the  critical 
damping  ratio  <£,  =  C,  /2nr\CJol  can  be  obtained  as  a 

function  of  amplitude.  The  is  found  to  be  constant  for 

small  amplitudes.  Some  values  measured  at  different  d/Jt  are  • 
shown  in  Fig.  18  along  with  a  typical  record.  At  d /l  =  0.059  a 
second  mode  damping  raolo  is  measured  and  is  found  to  be  ^  =  0.018. 


When  the  wind  is  turned  on,  at  low  velocities,  the  cable 
bends  back  and  assumes  an  inclined  static  equilibrium  position. 

Any  disturbance  imparted  on  it  quickly  dies  down.  An  eye  estima¬ 
tion  of  the  cable  inclination  is  made  from  outside  the  wind  tunnel 
using  a  transparent  protractor  fixed  to  one  of  the  clamped  ends. 
This  is  shown  in  Fig.  20,  for  d/J,  =  0.059  .  For  other  sag  ratios 

the  same  type  of  variation  is  observed  with  little  scatter.  The 
solid  lines  give  the  theoretical  values  calculated  from  Eq.  5.11 
which  will  be  discussed  later.  The  wind  velocity  is  increased 
in  small  steps  and  beyond  a  certain  critical  velocity,  the  cable 
begins  to  oscillate  by  itself.  The  amplitude  increases  until  a 
steady  state  value  is  reached.  The  mid-span  amplitude  is  meas¬ 
ured  by  a  vertically  travelling  telescope  from  outside  the  tunnel 
by  focussing  near  the  midspan  and  freezing  the  motion  of  the 
cable  by  a  strobe  light  which  is  also  used  for  measuring  the  fre¬ 
quency.  The  velocity  is  increased  in  small  steps  to  about 
100  ft/sec  and  then  decreased.  The  frequency  and  midspan  ampli¬ 
tudes  are  estimated  at  each  velocity,  and  are  shown  plotted  in 
Fig.  22a, b. 

At  low  velocities,  the  cable  first  begins  to  oscillate  in 
second  bending  mode  with  the  cable  at  blown  back  position.  As 
velocity  is  increased,  it  begins  to  oscillate  in  first  bending 
mode  in  a  swinging  motion.  At  a  higher  velocity,  the  second  mode 
appears  again.  The  cable  oscillates  primarily  in  the  first  mode 
with  small  second  mode  superimposed,  giving  the  appearance  of  a 
travelling  wave  across  the  span.  A  brief  outline  of  the  response 
description  is  given  in  Table  *1 .  At  mode  transition  points  and 
when  both  the  modes  occur  together,  the  second  mode  has  approxi¬ 
mately  twice  the  frequency  of  the  first  mode. 

It  is  observed  that  the  frequencies  measured  with  velocity 
increasing  and  decreasing,  follow  the  same  path  but  the  amplitude 
measurements  show  a  certain  amount  of  hysteresis.  The  arrows  in 
Fig.  22a  show  the  path  with,  velocity  increasing.  The  first  mode 
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frequencies,  nondimens ioned  by  w  ^  and  the  amplitudes  nondimen- 
sioned  by  sag  d  are  plotted  in  Fig.  23.  For  d/£  =  0.080,  0.059, 
0.043,  0.059,  0.043,  0.029  the  nondimensioned  frequencies  seem  to 
collapse  into  a  single  line  and  increase  almost  linearly  with 
velocity.  As  for  the  amplitudes  there  is  no  general  pattern.  For 
d/H  =  0.059  the  amplitude  is  quite  large  and  cable  swings  in  an 
almost  180°  arc  with  a  near  horizontal  mean  position.  For  higher 
d/£  =  0.080  the  amplitudes  are  lower  but  still  quite  large.  With 
reappearance  of  small  second  mode  near  70  ft/sec,  the  amplitudes 
tend  to  go  down.  For  lower  d/s,  =  0.043  the  amplitude  response  is 
similar.  With  reappearance  of  strong  second  mode  near  90  ft/sec 
the  amplitudes  go  down  rapidly.  At  very  low  d/2-  =  0.020,  0.011 
the  amplitudes  are  quite  large  compared  to  the  sag  and  increase 
rapidly  with  velocity. 

These  are  two  possible  reasons  for  the  irregular  amplitude 
responses  at  different  sag  ratios.  As  seen  in  Section  4,  the 
aerodynamic  forces  are  quite  sensitive  to  small  imperfections  at 
the  edges  of  the  cable  and  since  no  specific  side  of  the  square 
was  chosen  to  face  the  wind,  two  different  sides  facing  the  wind 
at  two  different  sag  can  give  unexpected  different  results.  Sec¬ 
ondly,  since  the  cable  is  not  entirely  free  from  small  initial 
twist  at  the  clamped  ends,  it  can  affect  the  amplitudes  to  a  cer¬ 
tain  extent. 

6.3  Theory  and  Experiment  Comparison 

The  static  aerodynamic  force  data  measured  for  the  square 
section  cantilever  beam  and  described  in  Section  3  can  be  used  for 
the  cable  since  both  are  of  same  cross-section  shape  and  the  ex¬ 
periments  are  in  the  same  Reynolds  number  range.  For  numerical 

calculation  an  average  representation  of  and  Cc  has  been  used. 
CN  =  1.0  to/n  p  —  10.0  tcvr^  p 

Cc  =  \.A  (6.1) 
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For  different  values  of  and  Qq  in  the  range  obtained  in  the 
test,  the  solutions  given  by  Eq.  5.21  changes  slightly  but  the 
nature  of  the  solution  remains  unchanged. 

To  check  the  static  representation  of  the  cable  one  can 
compare  the  linear  first  natural  frequency  at  zero  velocity  given 
by  Eq.  5.22.  This  is  simply  the  frequency  of  small  horizontal 
oscillation  of  the  sagging  cable. 

(Expt . ) 


d/S, 

iool  (Theory) 

H 

O 

3 

0.080 

2.5 

2.5 

0.059 

2.86 

2.8 

0.043 

3.3 

3.3 

0.029 

3.9 

4.0 

0.020 

5.2 

4.6 

0.011 

6 . 1 

5.5 

They  are  found  to  be  in  good  agreement  except  at  small  sag  ratios. 

The  static  positions  at  different  velocities  are  computed 
from  Eq.  5.11  and  plotted  in  Fig.  20  for  d/2-  =  0.059.  The  top 
picture  shows  the  side  view  of  the  cable  position,  and  these  are 
similar  to  those  observed  during  the  experiment.  The  inclination 
of  the  cable  from  the  vertical  is  also  shown  in  the  lower  picture 
and  is  found  to  compare  well  with  the  experimental  observations. 
For  other  d/S,  ratios  the  results  are  found  to  be  similar. 

The  steady  state  amplitudes  are  calculated  from  Eqs.  5.21 
by  harmonic  balance  method,  assuming  no  second  bending  mode  is 
present.  The  results  for  d/2,  =  0.059,  0.0*13  are  shown  in  Figs. 
24a, b.  An  estimation  of  maximum  induced  angle  of  attack  from  the 
frequency  and  amplitude  shows  that  it  is  well  within  the  linear 
Since  the  structural  damping  is  not  truly  represented 


range. 


the  solutions  with 


by  the  critical  damping  ratios  ^ 

£  ^  =0  are  also  shown  in  Pip.;.  24a.  Exclusion  of  structural 

damping  does  not  change  the  amplitude  substantially  except  near 
Jie  critical  velocity.  The  frequencies  agree  quite  well  with 
the  experimental  values.  The  amplitudes  are  only  in  qualitative 
agreement  for  these  two  sag  ratios.  Theoretical  results  for 
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other  sag  ratios  give  frequencies  which  are  close  to  the  experi¬ 
mental  values  and  the  amplitudes  follow  the  same  theoretical  trend 
as  for  d/l  =  0.059  and  are  in  disagreement  with  the  experimental 
values  except  for  the  similar  trend  of  large  increases  in  amplitude 
near  the  critical  velocity.  From  the  experimental  results  it  is 
observed  that  the  aerodynamic  lift  forces  do  not  reach  the  values 
given  by  the  quasi-static  theory  due  to  the  lack  of  spanwise 
correlation. 
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SECTION  7 


CONCLUSIONS 


The  present  report  has  investigated  the  self-excitation 
characteristics  of  a  square  section  body  in  the  form  of  a  vertical 
cantilever  beam  and  a  sagging  suspended  cable  when  placed  in  an 
uniform  wind  stream.  The  self-excitation  is  due  to  negative  damp¬ 
ing  type  aerodynamic  forces  induced  by  transverse  motion  of  the 
body.  Unlike  the  vortex  excited  oscillation,  this  kind  of  self¬ 
excitation  can  take  place  over  a  wide  range  of  velocity,  inde¬ 
pendent  of  the  natural  frequency  of  the  body. 

Since  no  purely  theoretical  method  of  predicting  the  aero¬ 
dynamic  forces  on  a  bluff  body  with  turbulent  wake  exists, 
experimental  quasi-static  representation  of  aerodynamic  forces 
was  used  to  determine  the  responses  theoretically.  They  were 
compared  with  experiment  to  examine  their  applicability  for  large 
three  dimensional  displacements. 

Experiments  were  conducted  with  a  square-section  flexible 
cantilever  beam,  free  to  deflect  in  any  direction  without  twist¬ 
ing,  and  placed  at  various  initial  angles  of  attack.  The  static 
aerodynamic  forces  wore  measured  for  use  in  the  quasi-static 
theory.  The  dynamic  responses  and  transient  motions  were  recorded. 
The  steady  state  amplitude  in  general  increases  monotonically  with 
the  wind  velocity  but  in  certain  cases  the  amplitude  drops  off 
over  a  range  of  velocity  (Fig.  11a).  In  these  cases  the  response 
was  found  to  be  very  sensitive  to  small  changes  in  initial  angle 
of  attack  aQ  (Fig.  12).  The  linear  equations  give  unbounded 
solution  beyond  a  certain  critical  velocity  and  the  amplitude  of 
oscillation  is  limited  by  nonlinear  aerodynamic  forces.  The 
theoretical  solution  for  the  steady  state  amplitude  agreed  only 
qualitatively  with  experiment,  generally  overestimating  the  values 
and  failing  to  predict  the  peculiar  amplitude  drop-off  behavior 
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(Pigs.  11a, b).  When  the  beam  was  at  an  initial  angle  of  attack  a0 
the  agreement  was  again  qualitative  (Figs.  12b,c,g,h).  The 
discrepancies  were  probably  due  to  imperfections  of  the  beam, 
unsteady  aerodynamic  forces  and  Reynolds  number  effects. 

Experiments  were  also  conducted  with  a  square  section, 
sagging  suspended  cable.  Large  self-excited  oscillations  were 
observed  and  the  amplitude  and  frequency  of  steady  state  oscilla¬ 
tions  were  recorded  (Fig.  22).  The  frequency  response  showed  a 
regular  trend  but  the  steady  state  amplitude  values  varied  widely 
with  different  sag  of  the  cable.  In  the  equations  of  motion,  the 
tension  in  the  cable  as  well  as  aerodynamic  forces  contribute  to 
the  nonlinearity  and  can  affect  the  steady  state  amplitude.  In 
the  first  mode  oscillations,  the  motion  remained  in  the  linear 
range  of  aerodynamic  forces  and  the  amplitude  was  limited  by  the 
tension  forces  together  with  the  damping  forces  due  to  the  drag. 
Reasonably  good  agreement  between  theory  and  experiment  w as  obtain¬ 
ed  for  the  frequencies,  but  the  steady  state  amplitudes  were  in 
qualitative  agreement  only  for  two  particular  sags  (Figs.  2^la,b). 
Again,  these  discrepancies  were  probably  due  to  imperfections  of 
the  cable,  unsteady  aerodynamic  forces  and  Reynolds  number  effects. 

Thus,  the  aerodynamic  forces  acting  on  an  oscillating  bluff 
body  with  sharp  edges  were  found  to  be  rather  complex  and  quite 
sensitive  to  small  imperfections  and  sharpness  of  the  edges, 
giving  rise  to  unexpected  responses.  The  actual  motion  induced 
forces  do  not  reach  the  values  given  by  quasi-static  theory  due 
to  imperfections,  lack  of  spanwise  correlation  and  unsteady  aero¬ 
dynamic  effects. 

The  results  obtained  here  may  be  of  interest  in  understand¬ 
ing  the  wind  induced  oscillations  of  tall  towers,  buildings,  launch 
vehicles  on  the  ground  and  power  lines.  To  understand  thoroughly 
the  unsteady  aerodynamic  forces  acting  on  bluff  bodies  with  sharp 
edges,  more  dynamic  tests  need  to  be  conducted. 
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APPENDIX  A 


CANTILEVER  MODE  SHAPE 


The  first  bending  mode  shape  of  the  uniform  cantilever 
beam  f^(5)  of  Eqs.  3.2  and  3.19  is  chosen  so  as  to  satisfy 
the  free  vibration  equation  and  all  the  boundary  conditions 
(Ref.  30)  .  ' 


|n^)—  Cosh  -  Cos  ^  -  o(r\  (  Sinhf*  ^  -  Sm  ) 

For  first  mode  <5,  =  1-8751  }  =  0.73-41  Sc  fi0)~2 
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A|=J'U(?)]%  =' 
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APPENDIX  B 


HARMONIC  BALANCE  COEFFICIENTS 
(FOR  CANTILEVER  BEAM) 

The  functions  f5,  (cu,  a*,  b2 )  e\c.  are  obtained 
from  the  integrals 

I  (a,fb„Qi4ba)  =  j\  2  Z  +  a,  Gw'c) 

-jr  m=3  v\*o  LL 

Vi 

x  U,  SmiT  4  ax  Gw't)  SintT  <k 

*  ~ 

fc,(a,,b„  a,A)  =  4j  2  £  %;„(-b,WC  +  a,G„r)w 

-or  tn=3  n=o  ^ 

M 

*  £  — 1>2  Sw  tr  4  a2G/it )  QonX 


and  similarly  for  Is 2  C  aL|  >  b,  »  a.?  ,  b2  )  and 

fC2  (a,,  b,  .  a*  ,  b2 )  by  replacing  ZmM  by  Ymr, 

If  b,  =  O  it  is  possible  to  get  a  general  ex¬ 

pression  for  these  integrals  as  given  below. 
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Case  II 


m  even,  n  odd 
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Hmny>  =  ^  [K(<-2)(K-4)  —  3,l][(f-Q(r-3) 

[(rfl+tt  +t)  (m-e-'n  —  l)  (w\+y\ -3  )  •  •• 


3.1] 
4.2  3 


He-ye  k  =  Vn  +  n-V~ 

If  any  factorial  or  first  bracketed  term  are  negative,  they  are 
set  equal  to  unity.  Also  01  =  J  . 

\  sz  (  >  0  j  0  b2)  ,  fez  CQ-i,  °j  a2  ^  )  are  given  by 

the  same  expression  after  replacing  2mn.  by  Ymn 
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APPENDIX  C 


FORCE  COEFFICIENTS 

The  quantities  Z*.  Yo#Y,/"Y5  used  in 

Kq.  3.46  are  defined  as  below  for  |  tamp  tom  |  ^  Q  [ 

Let  Pn  =  (  P«  —  Gy\  tan  ou )  Qj/3^oU 
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are  same  as  above  if  Po,^***^  are  replaced 
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APPENDIX  D 


CABLE  MODE  SHAPE 


The  mode  shape  {,(£)  y  f2  (?)  used  in  Eq.  (5.1*0 
are  the  first  and  second  bending  mode  shapes  of  a  freely 
vibrating  clamped-clamped  beam  (Ref.  30). 
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The  Amn  }  bai  ,  t>,  ,  t>2  /  Pr  are  mode  shape 

Integrals  as  defined  below. 
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5  mn  =  f  ^  frv  (?)  Je 

0°  d  f  d  £ 


68 


Dn  =  l'  l  (5)  ?„(?) 

-  &az/  Bu 


It  may  be  noted  that  =  o  has  been  used  in  the  formu- 


lation. 

&u  j  ,  etc. 

are  tabulated 

in  Ref.  30 

En, 

=  4.877 

D.  = 

138.46 

B22 

~  46.05 

D2  « 

3803.  5 

The  Anin  are  evaluated  by  numerical  integration  by  using 
Waddel's  rule. 


A20  =  °*^6 


1.0 


AiJ0  =  0,734  k22  =  1,:L1 

'  =  1.689  a'  =  0.44 

A  04 


Cn 


APPENDIX  E 


HARMONIC  BALANCE  EQUATIONS 

For  the  harmonic  balance  solution  of  Eq.  5.21  ,  it  is 
convenient  to  make  a  change  of  variable 

b  =  i  +  L 

Eq.  5.21  reduces  to  simply 

k  +  +  <X't!m(K2+  K2)  -  o  »>=  t,2 

where  '/  =  T\ 1  -  <X,  (  ) 

ti-  5, 

&  - 

For  convenience  all  ~  sign  are  replaced  by  a  first  subscript  2  . 
Thus 

b  -  a  =  b|0  +  +  O-n^yvva  CoC  +  b| z0oo  20>T  +  ai2Sin20/t 

ri  /f(0 

b-a  -  k  +  b2lGoU)t  +■  Qz,S'mO)t  +  b2ZGs»2Wt  4- azz5in  ZWC 
^  **2  0 

replaces  Eq.  5.23  . 

The  algebraic  equations  to  be  solved  are 
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Constants: 


£ 

“■'2  [CA.+  Q„,A„,+  *>„,>„*  + 

+  b*,»y  “  =  ° 

Cos  cat: 

-  2 

Y'-1 

■f2(bm,&ni+  Ow,  A„a  bw2&ni  +  QWJ.An,  4  bM1P>n3  +  Oms,An3)] 
Sin  cat: 

(y-  u2)  ami  -  <*) ^  bm,  +  tf,  £  {(Om  +  amtBno) 

«»l  ^  / 

+  2(bm,A«2+  a»*6m  -Qm&nZ  -bmaAnj  +  bm2An3  -  aWiP>*s)J 

Cos  2a?t: 

(V-^bvnj.  +  “%,<»,»*  +  «i2{(kmoB«+  bwi5n<,) 

+  z(bm|Bn|  -  am,  Arx  I  +  bmi5na+aTrMAYl2>+  b^  Bn^  +  Qyyi2L4y,4)J 

Sin  2c^t: 


('y-a^o^  -  cu|mbmi  + 

(*i  /L  Ani+  Q*i2.6>n0) 

n»i  1  J 

^  2  (  Am  +  Om 

^ni  +  ^mi  An*,  -o^b^  4-  bm2A«-q  ~  °vnz 

Wl  =  1,2 

where 

Co 

—  b  mo  •+■ 

B  ^0 

-  k*  Z  -V 
b  no  * 

2  C  bm  +  a*t  +  bni  -V 

&r\i 

=*  2b*0bn, 

+  bni  b^  2  +  °vu  an  -z 

A  ni 

=  2b*M6am 

■+■  bm  Qt)x  ~  brix^ni 

7' 

Bn-j.  — 

2  b  no  byv2. 

4 

H 

bru  - 

An  2.  — 

* 

2  bno  0  v\q. 

4 

t>  r\\ 

Oy\  i 

= 

bm  b  n  2. 

— 

Ovu 

0v»2 

Ar^j  = 

bm  Ona 

4- 

On> 

bni 

&r\  4  — 

i  (  b„\  - 

Qv\a  ) 

>  .  A  ^  — 

bm  0^2 

Newton  Raphson’s  iteration  scheme  for  solving  a  set 
of  algebraic  equations. 


f)  =  O  =  ■ 

r  ' 

l 

■  +  UlJ 

4X 

WHERE  [j]  = 

■  -I  YH-I 

n. 

n 

One  assumes  a  trial  solution 
solves  the  equations 


and  at  each  step 


[Jj. 


r 


1 

AX 

=  —  ■ 

i 

i 

X 

—  , 

| 

j* 

+ ! 

J 

r  \ 

AX 

. 

^  n  ! 

•  - 

n+i  1 

L  J 

n  1 

•.  - 

Convergence  criteria  A  X  |  <,.10^ 


Xj  is  of  order  1. 


1? 
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TABLE  1 


CANTILEVER  BEAM  PARAMETER 


Length 

1 

13-1  inches 

Side  of  square 

b 

0.09  inches 

Mass/length 

m  = 

0.01233  lb/ft 

Young's  modulus 

E 

30  x  10^  psi  (steel) 

El/mgl^  = 

9.58  (From  u,  ) 

= 

9.7^  (I  for  round  ro 

First  bending  mode  natural  frequency 


CPS  (Theory  9.5  CPS  with  I  for  round  rod) 


TABLE  2 


POLYNOMIAL  COEFFICIENTS 

The  following  table  gives  the  coefficients  Pn>  Qn  of 
8th  degree  polynomial  fit  of  CN,  Cc  vs.  TanS  curve  of  Fig.  7, 


8  respectively. 

CN  vs.  TanB 

Position  A 

A 

D 

Velocity 

37F/S 

65F/S 

37F/S 

p 

0.025 

0.025 

-.007 

o 

h 

0.88 

0.62 

1.35 

P2 

0.25 

0.07 

-0.7 

P3 

4 .  37 

7.67 

-14.7 

pt 

-10.0 

-3.9 

-4.3 

P5 

-151.8 

-181.6 

-18 . 3 

p6 

117.2 

52.0 

83.9 

e, 

1 

454.7 

526.2 

172.8 

P3 

-340.0 

-144.5 

-263.8 

CN  =  Z 
i\*o 

7« 


|fa*p  \  4 


TABLE  2  (Continued) 


Cc  vs.  Tan3 


Position 

A 

A 

D 

Velocity 

37F/S 

65F/S 

37F/S 

Qo 

1.2*1 

1.54 

1.07 

*1 

-0.3 

-0.47 

0.2 

Q2 

-9.8 

-14.9 

-3.1 

Q3 

7.9 

9.1 

1.1 

Q4 

200.1 

255.1 

115.5 

% 

-66.7 

-60.4 

-12.4 

Q6 

■1208.  4 

-1471.1 

-789.2 

Q? 

159  •  5 

123.1 

24.9 

Qg 

243*1.2 

2890.0 

1699.0 

°C=  L  ^  I  Up  U  CMS 


TABLE  3 


CABLE  PARAMETER 


Span 

Square  section  side 
Weight/length 
Bending  rigidity 
Tensile  stiffness 
E  I/mgl^ 

AE/rngl 


1  =  3.02  ft 
b  =  0.09  inch 
mg  =  0.00675  lb/ft 
El  =  0.0019  lb-ft2 
AE  =  7090  lbs 
=  0.01 

=  0.3^78  x  106 
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TABLE  4 


CABLE  RESPONSE  DESCRIPTION 

U  =  Free  stream  velocity  in  FT/SEC 

=  1st  mode  oscillation  frequency  in  CY/SEC 
w2  =  2nd  mode  oscillation  frequency  in  CY/SEC 
6  =  Inclination  in  degrees  (see  Fig.  20) 

W^/d  =  Midspan  single  amplitude/span  (see  Fig.  24) 
IBM,  2BM  -  1st  and  2nd  bending  mode 


a/i  =  o.oii 


u 

U_1 

e 

Vd 

Remarks 

0 

5.5 

- 

0° 

- 

- 

20 

- 

16.3 

o 

o 

m 

- 

Small  2BM 

26 

- 

17.5 

35° 

- 

ir 

29 

9.5 

19.0 

40° 

0.4 

IBM,  Small  2BM 

39 

13.2 

- 

■"65° 

0.95 

IBM 

51 

14.8 

- 

O 

o 

CO 

f 

1.02 

IBM,  Small  2BM 

71 

18.6 

- 

- 

1.6 

it 

61 

16.3 

- 

- 

1.3 

ii 

34 

11.0 

- 

'60° 

0.7 

n 

26 

_ 

18.0 

o 

o 

Small  2BM 

“  '•»-***•  ;ev.y,_»  w-msk  -  s 


TABLE  4  (Continued) 


d/t  =0.02 


u 

^1 

^2 

0 

Vd 

Remarks 

0 

4.6 

- 

0° 

- 

- 

29 

- 

13.8 

- 

Small  2BM 

36 

7.8 

lr.  3 

65° 

0.16 

Small  IBM  &  2BM 

39 

8.2 

17.0 

70° 

0.2 

it 

47 

10.3 

- 

78° 

0.28 

IBM,  Small  2BM 

60 

13-5 

- 

85° 

0.57 

it 

73 

16.0 

- 

- 

0.74 

ti 

81 

17.8 

- 

- 

0.85 

u 

59 

13.3 

- 

O 

o 

CO 

t 

0.58 

it 

45 

10.0 

20.0 

"75° 

0.21 

m 

27 

— 

13.0 

"50° 

.  , 

u 

8? 


TABLE  4  (Continued) 


0.029 


u 

fl 

_2 

6 

0 

4.0 

- 

0° 

32 

- 

- 

— 4 

O 

O 

45 

- 

- 

0 

85 

50 

8.3 

- 

- 

68 

11.0 

- 

- 

78 

12.5 

- 

- 

88 

14.2 

28.4 

- 

99 

15.7 

31.3 

- 

75 

12.2 

- 

- 

56 

9.2 

- 

~ 

50 

8.3 

- 

- 

42 

W^/d  Remarks 

Small  2BM 
Small  IBM  &  2BM 

0.2 3  IBM 

0.25  IBM,  Small  2BM 

0.12  IBM,  2BM 

0.06  2BM,  Small  IBM 

—  tf 

0.14  IBM,  Small  2BM 
0.24  » 

0.09  Small  IBM 

Small  2BM 


36 


78° 


Quiet 


TABLE  4  (Continued) 


0.043 


u 

^2 

6 

Vd 

Remarks 

0 

3.3 

- 

0° 

- 

- 

28 

- 

8.8 

57° 

- 

Small  2BM 

39 

- 

11.0 

72° 

- 

ir 

55 

7.5 

- 

O 

oo 

l 

Small  IBM 

58 

7.8 

- 

- 

IBM 

65 

8.6 

- 

- 

'  • .  5 

IBM,  Small 

76 

10.0 

- 

- 

0.6 

it 

88 

11.3 

- 

~ 

'>.65 

IBM  &  2BM 

100 

13.8 

- 

- 

0.43 

ii 

94 

12.3 

- 

- 

0.53 

it 

82 

10.5 

- 

- 

0.64 

it 

69 

9.2 

- 

“ 

0.6 

IBM,  Small 

54 

7.4 

- 

- 

0.49 

ti 

47 

— 

13.0 

o 

o 

CO 

1 

Small  2BM 

TABLE  4  ( Coiit inued) 


d/l  ■  0.059 


u 

c 

0 

Vd 

Remarks 

0 

2.8 

- 

0° 

- 

- 

29 

- 

- 

45° 

- 

Quiet 

41 

- 

- 

o 

o 

- 

Quiet 

61 

- 

13.8 

- 

- 

Small  2BM 

62 

6.9 

- 

- 

0.52 

IBM,  Small 

73 

7.7 

- 

- 

0 .66 

tf 

83 

9.1 

- 

- 

0.74 

It 

90 

9.8 

- 

- 

0.74 

If 

73 

7.9 

- 

- 

0.74 

II 

61 

6.7 

- 

- 

0.68 

II 

49 

5.5 

- 

- 

0.44 

It 

46 

- 

10.5 

~75° 

- 

Small  2BM 

36 

9.3 

~60° 

. 

it 

Or 


TABLE  4  (Continued) 


d/a  =  0.08 


u 

^1 

0 

2.5 

- 

38 

- 

7.7 

47 

WE 

9.2 

50 

4.7 

- 

63 

5.8 

- 

70 

6.5 

- 

83 

7.8 

99 

9.3 

93 

8.8 

- 

76 

7.1 

- 

54 

4.9 

«n 

47 

- 

9.2 

24 

r. 

«• 

0 

Vd 

Remarks 

0° 

~ 

- 

62° 

- 

Small  2BM 

75° 

- 

it 

o 

CO 

P"— 

0.36 

IBM 

o 

o 

CO 

0.51 

tt 

85° 

0.52 

IBM,  Small  2BM 

- 

0.48 

n 

«c. 

0.31 

tt 

- 

0.33 

tt 

- 

0.44 

tt 

~ 

0.34 

tt 

o 

0 

- 

Small  2BM 

30° 

t» 

86 
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66  FT/SEC 


FIG.  2a  SCHEMATIC  DIAGRAM 


FIG,  2b  COORDINATE  SYSTEM 


FIG.  4.  VELOCITY  AT  SPECIMEN  POSITION  VS.  REFERENCE  VELOCITY 


IG.  5a  CT  VARIATION  WITH  STATIC  ANGLE  OP  ATTACK  a  (POSITION  A) 


VELOCITY 


VELOCITY 


VARIATION  WITH  STATIC  ANGLE  OP  ATTACK  a  (POSITION  D) 


FIG.  7b  C„  VS.  tan 6  (POSITION  D) 


VS.  tans  (POSITION 


VS.  tans  (POSITION  D) 
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FIG.  9  STRAIN  GAGE  CALIBRATION  CHART 


FIG.  lib  STEADY  STATE  AMPLITUDE  RESPONSE  (THEORY  &  EXPT 


WT/b  w-r/b 


STEADY  STATE  AMPLITUDE  AND  INCLINATION  AT  INITIAL  ANGLE  cx0  (POSITION  A) 


FIG.  12c, d  STEADY  STATE  AMP.  ’TUDE  AND  INCLINATION  AT  INITIAL  ANGLE  a0  (POSITION  A) 


STEADY  STATE  AMPLITUDE  AND  INCLINATION  AT  INITIAL  ANGLE  aG  (POSITION  A) 


12g,h  STEADY  STATE  AMPLITUDE  AND  INCLINATION  AT  INITIAL  ANGLE  a0  (POSITION  D) 
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Wr/b 


-10 


5 


1 


10  SEC  - H 


U  =  I7.3FT/SEC 
w  =  3.4  CPS 

:?< .  \  o‘t'.  M  No  1. 1 
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Wr/b 

.10. 


.  Wt/o 


10 

5 

0 


■wei: 


10  SEC  — I 


U  =  37  FT/ SEC 
o  00=  9.4  CP6 


FIG.  13  SAMPLE  TRANSIENT  RECORDS  (POSITION  A,  a0  =  0) 
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FIG.  15  LINEAR  THEORY  -  EQUIVALENT  -  DAMPING  RATIO 


!  IG .  16  (  CONT'D  )  GRAPHICAL  REPRESENTATION  OF  ST  RADY  STATE 
SOLUTION  AND  STABILITY  FOR  1  D.O.F.  MOTION 
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